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Abstract. Fermi Normal Coordinates (FNC) are a useful frame for isolating the locally 
observable, physical effects of a long-wavelength spacetime perturbation. Their cosmological 
application, however, is hampered by the fact that they are only valid on scales much smaller 
than the horizon. We introduce a generalization that we call Conformal Fermi Coordinates 
(CFC). CFC preserve all the advantages of FNC, but in addition are valid outside the horizon. 
They allow us to calculate the coupling of long- and short-wavelength modes on all scales 
larger than the sound horizon of the cosmological fluid, starting from the epoch of inflation 
until today, by removing the complications of the second order Einstein equations to a large 
extent, and eliminating all gauge ambiguities. As an application, we present a calculation of 
the effect of long-wavelength tensor modes on small scale density fluctuations. We recover 
previous results, but clarify the physical content of the individual contributions in terms of 
locally measurable effects and “projection” terms. 
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1 Introduction 


When performing a calculation one may choose any set of coordinates. However, choose the 
“wrong” coordinates, and you will be sorry. 1 This is especially true in general relativity 
and its application to cosmological perturbation theory. A first and most obvious issue is 
that the “wrong” coordinates might substantially increase the computational complexity of 
the problem. A second and maybe more subtle issue is the risk of incorrectly interpreting 
a legitimate solution of the equations of motion in terms of measurements because the co¬ 
ordinates in which that solution has been derived are different from the natural coordinates 
of the observer. One often loosely describes this situation as gauge artifacts. Of course, the 
computation of a physical observable, defined in an explicit operational sense, must give the 
same result in any set of coordinates. However, even the simplest operational measurement 
might become very subtle to define when using the wrong set of coordinates. In the context 
of cosmology, Maldacena’s consistency relation [1, 2] provides a sharp example of how subtle 
gauge artifacts can be. The consistency relation states that, for single clock inflation, pri¬ 
mordial non-Gaussianity is of the local form with oc (n 3 — 1), where n s is the spectral tilt 
of the primordial scalar power spectrum. One might naively think that the locally measured 
short-scale scalar power spectrum in the presence of a long scalar mode (say of () is different 
from what one would measure locally in the absence of the long mode. One might even be 
tempted to speculate that if (n s — 1) were allowed to be somewhat larger than in standard 
slow-roll inflation, one could have a significant locally observable effect. This interpretation 
and speculation are both incorrect because the consistency condition is derived using the 
comoving coordinates of an unperturbed universe, while the ruler of a local observer would 
instead change in the presence of the long mode in precisely such a way that she would not 
see any difference from an unperturbed universe (see [11] for a detailed discussion). In fact, 
the consistency relations discussed widely in the literature [1-10] encode the gauge freedom 
or diffeomorphism invariance remaining in the description of single-field slow roll inflation. 
Single-clock inflation does lead to some observable non-Gaussianity of the local shape in say 
the Cosmic Microwave Background (CMB) and Large Scale Structure. However, this is due 
to the effect of large-scale perturbations (within the current horizon) on the propagation of 
photons from the source to the Earth (“projection effects”) which are entirely independent 
of interactions during inflation and do not scale as n s - 1. 

Choosing the natural set of coordinates of a local observer is very convenient because 
they are directly related to local measurements, hence explicitly eliminating any gauge arti¬ 
facts. In most applications to cosmology, the local observer is in an inertial frame, free-falling 
in the local gravitational potential. Such an inertial observer can describe the neighborhood 
around her worldline as flat spacetime with corrections that grow with the square of the 
distance from the observer’s geodesic times second derivatives of the spacetime metric as 
described by the Riemann tensor. This holds in an arbitrary spacetime, and the set of co¬ 
ordinates is known as Fermi Normal Coordinates (FNC) [12] and it is implicitly used by all 
laboratory experiments for which the gravitational field of the Earth is negligible, such as for 
example at the Large Hadron Collider. A major drawback that hampers the usage of FNC 
in cosmology is the fact that the Friedmann-Lemaitre-Robertson-Walker (FLRW) metric has 
second derivatives of order H 2 , where H is the Hubble parameter. In such a spacetime, FNC 
cover at most a patch of size H which is too small to describe horizon-scale or super-horizon 


1 This quote is allegedly due to S. Weinberg. 
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modes k < aH. Hence one cannot use FNC to connect primordial perturbations generated 
during the early universe to the late time evolution leading up to the time of observation. 

In this paper we formalize the construction of a generalization of FNC, which we call 
Conformal Fermi Coordinates (CFC). 2 These are the coordinates of a local observer that 
describes the spacetime in a neighborhood of her worldline as an FLRW spacetime. The 
corrections from the unperturbed FLRW spacetime again grow quadratically in the distance 
from the worldline. However, depending on the structure of the spacetime, the corrections 
can stay small on scales much larger than the Hubble horizon. Thus, these coordinates 
share all the advantages of FNC but are valid on super-horizon scales. Indeed, just as 
particle physicists implicitly use FNC, as cosmologists we implicitly use CFC in describing 
our background universe , allowing us to ignore deviations from the FLRW spacetime that are 
beyond our current horizon (even though they could be very large). 

The perks of using CFC in cosmology are best described for a toy spacetime that consists 
of just one long- and one short-wavelength perturbation on top of a homogenous and isotropic 
background. These perturbations might be scalars, vectors, tensors or a combination of them. 
The following considerations then apply: 

• The CFC defined for an observer that is inertial with respect to the long mode can 
be used to describe the dynamics of the short mode. At the observer’s location, it 
describes a local FLRW universe with tidal corrections, whether either or both the long 
and the short mode are inside or outside the Hubble horizon, as long as the long mode 
is outside the sound horizon of the fluid. For instance, in a universe dominated by dark 
matter, a spherical overdense region locally collapses in exactly the same way as an 
over-critical universe. Regardless of the size of the region compared to the horizon, the 
picture applies as long as the region is larger than the free-streaming distance of the 
dark matter. 

• As it is the case for FNC, results computed in CFC have a very transparent physical 
interpretation: They correspond precisely to what a local observer would measure, as 
she would find herself living in a FLRW spacetime whose Hubble parameter matches 
the locally-measured rate of expansion and whose spatial curvature matches the local 
spatial geometry. In this sense, CFC are the complement of consistency relations, in 
that they remove the diffeomorphim invariance to isolate all locally observable effects. 

• Using CFC and given some primordial perturbations, one can organize the computa¬ 
tion of observables in a convenient and physically transparent way (see e.g. [11]) by 
transforming to CFC at some point during inflation and then following the evolution 
of small scales in this frame. In a further step, one computes the so-called projection 
effects to relate local physics to what a distant observer sees. This is done perhaps 
most simply by transforming to FNC first (which takes into account the local CFC 
scale factor) and then using the ruler perturbations of [13] (Sec. 6). 

• The computation of non-linear effects that couple long and short modes is substantially 
simplified. This in particular applies to the case where the short-wavelength modes are 
far inside the horizon (which is the case typically of practical interest). We show here 
that it is sufficient to include the much simpler non-linearities in the fluid equations 

J They were first introduced in a less formal way in [11]. We make the connection to their construction in 
Sec. 4. 


- 3 - 



while keeping the Einstein equations to linear order in the small-scale modes. This had 
already been used in [16]. We also show how the CFC calculation can be extended 
to the case in which the short mode is super-horizon. In this regime one in general 
needs to include certain non-linear terms in the Einstein equations, and some subtleties 
arise in dealing with spatial derivative operators. In this regime, which is typically 
not phenomenologically interesting, the calculation is not any simpler than a global 
calculation, but it does retain the advantage of the more transparent interpretation. 

The remainder of the paper is organized as follows. In Sec. 2 we set the stage with 
a brief review of FNC. We then describe the construction of the CFC in Sec. 3. In Sec. 
4 we discuss how to choose a local scale factor appearing in CFC for a perturbed FLRW 
universe, and give formulae for the conversion to CFC. We also discuss the subtle differences 
to [11]. In Sec. 5 we derive the relevant terms in the Einstein and fluid equations needed 
to capture the leading-order couplings between long and short modes. This provides one 
of the main motivations for using CFC. We discuss in detail the relation to observations 
made by a distant observer in Sec. 6. Sec. 7 discusses application of the CFC to a long 
tensor mode affecting short scalar modes, extending the results of [16] to the case in which 
the short mode is super-horizon, and clarifying the interpretation of their result. We finally 
conclude in Sec. 8. Useful mathematical results are collected in the appendices. A more 
detailed treatment of long-wavelength scalar perturbations will be presented in a companion 
paper [29], 

2 Recap of Fermi Normal Coordinates (FNC) 

We first briefly review the basic concept of FNC and its construction before generalizing to 
CFC in the context of cosmology. 

Consider a free-falling observer, whose trajectory is a time-like geodesic, which we call 
the central geodesic. One can choose that it defines the spatial origin of the coordinate 
system at all times, and further that the tangent vector to the geodesic defines the time 
direction. The most obvious choice is simply to take the proper time tp along the central 
geodesic as time coordinate. We then choose an orthonormal set of spatial basis vectors, 
completing the tetrad (ea)^, such that at some point on the geodesic g^ u = r] FU . Parallel 
transport of these basis vectors ensures that this is valid for any point along the geodesic. 
Further, the fact that the observer is free-falling allows us to choose coordinates so that for 
any point along the central geodesic {dpg^p = 0. Thus, all Christoffel connections vanish, 
(r(L)_F- For points off the central geodesic, the FNC metric differs from the flat one only at 
quadratic order in the spatial FNC coordinate x l F . In other words, the FNC are the natural 
coordinates a free-falling observer would use to describe local measurements. The leading 
gravitational effect she would observe is a tidal field. 

Since the FNC is only necessarily flat in the vicinity of the central geodesic, the param¬ 
eterization of the spatial coordinates is not unique (see e.g. [11]). However, a conceptually 
simple geometrical construction of FNC satisfying all of the aforementioned requirements 
exists [12], The slice of simultaneity containing a given point P on the central geodesic is 
defined as the surface spanned by spatial geodesics that radiate outward from P. The mag¬ 
nitude of the FNC position x l F can be fixed by the proper length from P along the spatial 
geodesic. To be more specific, the spatial geodesic connecting P (having FNC coordinates 
x F (P) = {^,0}) and another point Q (having FNC coordinates x F (Q) = {rp,x F }) on that 
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slice can be parameterised by an affine parameter A, where A = 0 at P and A = 1 at Q. The 
direction of the geodesic is fixed by the initial condition 

= x l F (eiY P , ( 2 . 1 ) 

A=o 

where x F is taken to be the FNC position for Q. This ensures that the proper distance 
squared is given by 5ijx' l F x :) F in the vicinity of P, since gpv(ei)p(ej) u p = Sij. Under this 
construction, the coordinate transformation from some coordinate system x^ to the FNC 
coordinate x F can be computed order-by-order in x l F . One repeatedly uses the geodesic 
equation to determine the coefficients of the series expansion x^(A) in A, which are completely 
determined by derivatives of the metric evaluated at P. Truncation at third order is sufficient 
to determine the transformation Jacobian matrix dx^/dxp to second order in x F and the 
tidal corrections in the FNC metric. 

The FNC metric can then be found order-by-order in x l F using 3 

F dx a dx@ . 

9pv\ x F ) — q x h q xL , 9ap{%) ■ ( 2 - 2 ) 

F F 

Here and throughout, a sub- or superscript F denotes quantities in the FNC or CFC frame. 
One can then verify that this construction successfully realizes that the FNC metric is flat 
up to second order in the spatial deviation from the central geodesic, 


~d\ 


9^{xf) = + 0[R^(P F ) 2 ] , 


(2.3) 


where the magnitude of the corrections is given by the Riemann tensor (transformed to the 
FNC frame). As an example, we can consider a Friedmann-Robertson-Walker spacetime, 
given by 


ds 2 = o 2 (r) 


-dr 2 + ( 1 + ^A'x 2 


-2 


8ijdx l dx 3 


(2.4) 


where x 2 = 5ijX l x 3 and K is the curvature constant. This metric becomes in FNC (e.g., 

[17]) 4 


F 

9m — 

F 

9ij — 


- 1 + 



H2(tF) + YF) x2p 


F n 

9oi 0 j 

8ij , 


(2.5) 


where H = a~ 1 da/dt = H/a. Of course, one can derive a more general expression allowing 
for a perturbed FLRW metric. However, it is already clear from Eq. (2.5) that the FNC in 
the cosmological context are only valid on scales that are much smaller than the horizon, 
since we are expanding perturbatively in Hx l F \ if this quantity becomes order one, the per¬ 
turbative description of the FNC metric breaks down. We will now show how the conformal 
generalization of FNC can get around this limitation. 

3 Note that the coordinate transformation enters not only through dx a /dxp, but also through the explicit 
coordinate shift from the central geodesic, g a p(x) = (g a p) P + {g a p lP ) P fix' 1 + • ■ •, where starts at linear 
order in x'p. 

4 There is residual gauge freedom in the spatial metric components as discussed in [11, 29]. 
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3 Conformal Fermi Coordinates (CFC) 

We now describe the construction of the conformal generalization of FNC. Since we shall not 
discuss FNC and CFC at the same time in the following, the subscript F will be reserved for 
CFC hereafter. The CFC are constructed in the vicinity of a timelike central geodesic, just 
like the ordinary FNC. However, unlike the FNC, it does not restrict the local spacetime to 
be Minkowski, but allows for a homogeneous expansion over time; that is, the lowest order 
CFC metric is an FLRW spacetime Eq. (2.4). Thus, the CFC metric takes the following form 

&( x f) = a F( T F) + K^ Tf ’ x f)] » Ku = °[( x f) 2 }- C 3 - 1 ) 

Like in the case of FNC, corrections to the conformally flat part start at quadratic order 
in x l F . However, a few additional subtleties arise. First, the CFC time Tp should be some 
suitable conformal time rather than the observer’s proper time. Besides, a suitable local 
scale factor ap{jp) should be defined in a physical way. In particular, if g ^ describes an 
unperturbed FLRW metric (but given, for example, in some unusual coordinates), then the 
CFC construction should yield the metric in the canonical FLRW form, Eq. (2.4). These 
considerations motivate our generalization of the usual FNC construction presented in the 
previous section. 

3.1 Constructing CFC 

The geometrical relation between the CFC and the global coordinates is sketched in Fig. 1. 
Throughout this paper, “global coordinates” refers to some set of coordinates valid at least 
in the region surrounding the geodesic considered (for example in one of the standard gauges 
of cosmological perturbation theory). This name is chosen to distinguish it from the local 
construction of the CFC and FNC. 

First, we choose the same set of orthonomal tetrads (e(()p and (e^)p, i = 1,2,3 as in 
the construction of FNC. We also parameterize the observer’s geodesic, which is also the 
spatial origin of CFC, in terms of the proper time tp in the usual way. 

Next, consider a spacetime scalar cif(x), which we require to be positive at least in a 
finite region around the central geodesic. We can then define a “conformal proper time” Tp 
through 

drp = a~ F 1 {P{tp))dtp , (3.2) 

where P(t,F ) is the point along the central geodesic at proper time tp. This can be integrated 
to yield a unique relation Tp{tp) (up to an integration constant which can be absorbed into 
a redihnition of clf). We then choose Tp as our time coordinate. We will often write clf(tf ) 
instead of ap(P(Tp))-, however, one should keep in mind that clf and the Tp(tp) relation 
depend on the specific geodesic under consideration. 

Now we need to define the slices of constant Tp , which requires some care. They should 
not be defined by simply tracing out spatial geodesics orthogonal to the central geodesic. In 
fact, for a homogenous flat FLRW spacetime, parameterized by the conformal coordinates, 
straight lines on (3-dimensional) constant-conformal-time surfaces are not true geodesics, but 
are only geodesics with respect to the conformal metric, i.e. = a -2 #^. This suggests 
that surfaces of constant-rj? in CFC should be spanned by space-like conformal geodesics, 
namely geodesics with respect to the conformal metric 

g F „{x) = a F 2 (x)g lliy (x) . (3.3) 
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Note that for a perturbed FLRW metric = a 2 (j] ll v + h tw ). in general differs from 
rg w + lifw because the local scale factor a F is not necessarily the same as the global one 
a. Here and throughout, a tilde denotes quantities defined with respect to this conformal 
metric. Let us then summarize how to locate the point Q corresponding to CFC coordinates 
{t f ,x 1 f }: 


1. Find the point P on the central geodesic corresponding to the CFC time r F . This point 
thus has CFC coordinates {tf,0}. 

2. Let h{TF] cd; A) denote the family of conformal geodesics with respect to g liu , with the 
affine parameter at P given by A = 0 and the tangent vector at P given by cd( e i)p- 
That is, a 1 are constants specifying the inital direction of the geodesic while A measures 
the geodesic distance with respect to the conformal metric (up to a constant factor). 

3. The point Q is located on the conformal geodesic h(rp\aF(P)P l ] A) where A = ( 5ijX l F x 3 F ) 1 / 2 
and (3 l = x l F /8ijX l F x? F . This ensures that the proper distance squared from P to Q 

is a 2 F 5ijX l F x 3 F at lowest order, which is the desired relation in the metric Eq. (3.1). 5 * 


This uniquely specifies the CFC coordinates, which are guaranteed to be regular in a 
finite region around the central geodesic (note that a F > 0 is a necessary condition). Clearly, 
in order to properly specify CFC, it is not sufficient to have clf(t f) given only on the geodesic. 
Rather, we also need its spacetime derivatives away from the geodesic. This is because we 
need to integrate the equation for the conformal geodesics in order to reach a point away 
from the geodesic, which is given by 


d 2 x tl ~ u dx a dx 13 
__ = n 

d\ 2 a ? dX dX 


(3.4) 


Here the Christoffel symbols for the metric g^ u are replaced by the conformally transformed 
ones 


f U = ~ c “-f> ■ < 3 - 5 > 

The tensor of shift is given by [15] 

C^a /3 = So Vfi hi a F + S'g V a In a F - g a /3 9^ Va In a F . (3.6) 

In order to obtain the leading correction to the CFC metric, it is sufficient to specify the 
first and second derivatives of a F only. The first derivatives of a F are already constrained 
by Eq. (3.1). In order to obtain the CFC form, the gradient of a F along the central geodesic 
has to be along the time direction, i.e. 


In a f 


4 =° 


(In of)'of (e 0 ) M , 


(3.7) 


where a prime denotes a derivative with respect to t F - Note that the second derivatives 
V^VylnaFlcci^o are already fixed by the CFC construction, since in CFC a F = of(tf). In 
particular, they exactly match the expression for an unperturbed FLRW universe (App. A). 

5 Note that even for a curved FLRW universe, the proper distance is given by this relation to lowest order 

in xf- 
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Figure 1: Construction of the CFC. At point P, the observer’s geodesic G intersects a 
spatial hypersurface X, having constant conformal time r and scale factor a(r) in some 
global coordinate system. The spatial hyper-surface X^, having constant CFC conformal 
time Tp and CFC scale factor clf(tf), also intersects G at P, but does not coincide with 
X in general. Another point Q on X^? is connected to P by a conformal geodesic, and is 
parameterized by the same Tp but nonzero x l F . The CFC coordinates are valid within a 
tubelike region bounded by hyper-surface B surrounding G. 


With this construction, we can now derive the explicit transformation law from some 
global coordinate system to CFC. The conformal geodesic which satisfies Eq. (3.4) has a 
power expansion in the curve parameter A, 

OO 

x^X) = Y,< X "- ( 3 - 8 > 

n =0 

This curve connects point P with CFC coordinates x F (P) = {tf\0} with point Q which 
in CFC has the coordinates x F (Q) = {tf,x 1 f }. Since P is chosen as the spatial origin, we 
immediately have «q = x IJ/ (P). Rescaling A so that it runs from A = 0 at P to A = 1 at Q, 
the tangent vector at A = 0 is specified by x F through 

= a F (P){eiY p x l F , (3.9) 

A=0 

Higher-order coefficients otn can then be recursively computed using Eq. (3.4). This is shown 
up to third order in App. B. 


Oil = 


dx 11 

~d\ 
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3.2 CFC metric 


We now derive the explicit form that the metric attains in the CFC frame. For the ap¬ 
plications we will describe below and in [29], we are also interested in deriving the leading 
correction h^ w to the CFC frame metric. 

The CFC metric of course still obeys the usual tensorial transformation rule Eq. (2.2). 
This transformation law, evaluated at some point Q not necessarily on the central geodesic, 
can be recast into the following form: 


[(M-P)) 2 g F „\ 0 4) 


a(P) dx a 

a(P ) dx@ 

(a(Q)V r 

ap(P)dx F 

clf(P) dx u F 

u (p)j y 


9a/3 


(x A ), (3.10) 


so that it describes a transformation law for the conformal metric a~ 2 g fiu . The notation for 
the global scale factor is understood as a(P) = a(r(P)) and a(Q) = a(r(Q)). It is meant to 
be the background scale factor evaluated at the point of interest, and therefore only depends 
on the corresponding global time coordinate r. This however is different for P and Q, so 
that a(P) / a(Q): different points on a constant -tf surface (and hence corresponding to the 
same (1f{tf) in CFC) are not simultaneous in the global coordinates, as evident from Fig. 1. 
Therefore, while clf{Q ) = clf(P) because they refer to the same CFC time and hence the 
same CFC scale factor, a(Q) and a(P ) are numerically different. 

Eq. (3.10) can be straightforwardly evaluated using the coordinate transformation de¬ 
rived in App. B. However, this becomes quite lengthy and is not very illuminating. We instead 
present a simpler derivation through a generalization of the method adopted in Ref. [12]. The 
underlying idea is to project various geometric quantities into CFC, in which their compo¬ 
nents take simpler forms than in arbitrary global coordinates. 

In CFC, the temporal coordinate vector is (eo) = d/dx° F = a,F(P)(eo)p, and the spatial 
coordinate vectors are (e^) = d/dx l F = ap{P){ei)p. Hence the CFC metric right on the 
central geodesic must be conformally Minkowskian, 


9\ 


fiu\p 


(oX {e v f g a p J p = a z F (P) [(e F ) a ( e v ) p g af} J p = a z F {P) rj^, 


\P 


(3.11) 


as the tetrad vectors are orthonormal. Further, the spatial CFC coordinate lines are geodesics 
of the conformal metric, and are thus simply parametrized through 


x°p = Tp = const., x F = /3 l X , 


(3.12) 


where /3* are constants and A is the affine parameter. Since the tangent vector is (0,/T), the 
conformal geodesic equation Eq. (3.4) reduces in CFC to 


(f 


F\F 
V 


FF = 0, 


(3.13) 


which implies (f 77 )C | p = 0 because (3 1 is arbitrary. One can further show easily (App. C) that 
following Eq. (3.12), all Christoffel connection coefficients (r F )([ fl , computed with respect to 
g FV [Eq. (3.3)] and in CFC , vanish on the central geodesic, i.e. (r F )^|p = 0. It follows that 
all first-order derivatives of the conformal metric g in CFC vanish on the central geodesic, 


(■ d a g F v) F 


= 0. 
p 


(3.14) 


This proves the absence of 0[x^] terms in the CFC metric. 

















Furthermore, the 0[(xp) 2 ] terms in the CFC metric are found to be related to the 
conformal Riemann curvature tensor R , i.e. the Riemann tensor constructed for the metric 
Eq. (3.3), and evaluated on the central geodesic. Given Eq. (3.11) and Eq. (3.14), the entire 
derivation in Ref. [12], which is for an analogous relation between the FNC quadratic metric 
terms and the ordinary Riemann curvature tensor, can be borrowed over, once conformally- 
related quantities replace the ordinary ones everywhere (e.g. metric, geodesic equation, 
geodesic deviation equation, Christoffel symbols, Riemann curvature tensors, and so on). 
The end result for the quadratic corrections reads 


9oo( x f) 


9oi(x F ) 


afj{x F ) 


a 2 F (r F ) 



a F (r F ) 
a F (r F ) 


g -^0 kil 


$ij „ R 


X pX p 


ikjl 


k 1 

sy‘ r ' J rf» L 

X pX p 


(3.15) 

(3.16) 

(3.17) 


Here R F is the Riemann curvature tensor constructed with respect to g^ and transformed 
to the CFC frame. In terms of some global coordinates, 

Ra/3~/5 = (®«)^(®/l) Rpupa, (3.18) 


where Rpvpa is the Riemann tensor of the conformal metric with its components computed in 
the global coordinates. Note that the indices of the conformal Riemann tensor R are always 
to be raised and lowered with the conformal metric g 

Extending the analysis to higher orders, as detailed in App. C, the 0[(x' l F p\ corrections 
to the CFC conformal metric scale as 


(diRpjvk) p x' F x? F x k F , (3.19) 

i.e. they are suppressed by the spatial derivative of R F multiplied by x F . Further, the 
coefficient of the ©[(x^) 4 ] terms in general contains didjR F as well as ( Rp ) 2 . This illustrates 
that we are expanding in the spatial variation of the tidal force R F induced by the long 
mode, as well as its amplitude. 

3.3 Choosing the CFC scale factor 

So far we have not written down any equation that determines the spacetime scalar a F (x). 
Note however that since in CFC a F can only depend on x l F at second or higher order, and 
since we fixed the second order term by matching to FLRW solution, we have effectively 
reduced this freedom to a function a F (r F ) of time only. Thus, another way of saying this is 
that the aforementioned construction is invariant under the most general reparameterization 
of the conformal time, 


dr' 

Tp = t f (t f ), a F (r F ) = a F( T F)^T! ( 3 - 20 ) 

under which a F {r F )dr F = dt F is still the proper time interval along the central geodesic. The 
new x z F s then have to be re-defined following our procedure using conformal geodesics, which 
in general leads to a different slicing of simultaneity, and different CFC metric perturbations 
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Eqs. (3.15)-(3.17). Note in particular that we can simply set of = 1, in which case we 
recover the standard FNC frame described in Sec. 2. 

Let us investigate this issue a bit further. Consider a congruence of geodesics in the 
vicinity of the chosen central geodesic, and let tP = dx^/dtp denote the tangent vector 
to this congruence. Note that (eo)^ = EP defines the time component of our tetrad. U^ 
satisfies the geodesic equation U a \7 a U M = 0, and hence (CPV Q CP) = 0. Straightforward 
manipulation allows us to derive 

u a v a {v F u») = - (V M CP) (v q cp) - R^ a u p u a . (3.21) 

Because U^ is geodesic, the velocity shear tensor V ^U a is in the “velocity-orthogonal” 3- 
dimensional subspace, whose projector is given by 

Pt = K + U»U V , (3.22) 

noting that we use the mostly positive metric convention; explicitly, the geodesic equation 
implies that P^P^XPt/^ = \PtP. We now decompose the velocity shear tensor in the usual 
way as 

= + < + (3.23) 

where d = P^V^tP = V^tP is the velocity divergence, a is the trace-free symmetric velocity 
shear, while ui is the antisymmetric part, i.e. vorticity. That is, 6 a a p = V( ;i P a j — dP a p/3. 
Eq. (3.21) then becomes, using [PV Q = d/dtp and P a /gP a ^ = 3, and using the antisymmetry 
of Ul a 0 • 

= ~ <Ta ^ cra!j + 0Ja ^ uja ’ 3 ~ RlJ ' 0 iio Uplja ■ (3.24) 

This is the well-known Raychaudhuri equation, where typically the contraction of the Rie- 
mann tensor is replaced with the Ricci tensor. We will not do this because we prefer results 
in terms of the Riemann tensor. Note that Eq. (3.24) is a purely geometric result; we have 
not made use of Einstein’s equations. Working in CFC, we have CP = a i? 1 (l, 0,0,0), so that 
Eq. (3.24) becomes 

+ u ij u ' 3 ~ a F 2 ( RF Tou.o • ( 3 - 25 ) 

Using the transformation law of the Riemann tensor under a conformal rescaling of the 
metric, g FV —> g^ = adp 1 g liv . and some algebra we obtain 

*- F \ rF )X o = «f 2 (^ F )Vo - 3 (Hp + Hi ). (3.26) 

Here we have used ()' = ap(') when acting on a simple function of time along the geodesic. 
We can also use the fact that by way of the CFC construction, (P F )° 0 ^ 0 = 0 (this is because 
the tetrad is parallel-transported; see Sec. 6.1). Finally, we obtain the Raychaudhuri equation 
in CFC: 

A-d + ^d 2 = 3 {Hp + H 2 p) - a ij( T ij + t OijJi - a F 2 (R F ) j 0j0 . (3.27) 

throughout we are always raising and lowering indices with g^, not P flv . 
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Again, we have made no assumption about the spacetime or about clf so far. We only used 
the requirement that U^ is tangent to a congruence of geodesics. In other words, for any 
spacetime we can always choose coordinates where the Raychaudhuri equation takes this 
form. 

Let us drop a and uj for the moment (in the case of perturbed FLRW discussed in Sec. 4, 
they are at least linear in perturbations, so that their contribution to Eq. (3.27) is at least 
second order and hence will be dropped in our linear treatement there). We then have 

Atf + = 3(i ^ + H 2 p) _ a -\RFy 0jQ . (3.28) 

Recall that once we have picked our central geodesic, this equation holds at all points along 
that geodesic, and all quantities are just functions of time. We can now use our freedom 
to choose aF(tp) as we desire. For example, choosing clf = 1 reduces B T to R F , which 
then contains the entire source of local velocity divergence. Our goal is to absorb the leading 
contributions from the Hubble expansion in of however, and make R F “as small as possible”. 
Guided by this principle, we will define clf through 

H F = ^ := L (3.29) 

dtp 3 

throughout this paper and the companion paper [29]. This is a convenient choice, first of all 
because it reduces to the true scale factor in the case of an unperturbed FLRW spacetime, 
in which case R F then vanishes entirely. Second, id is then an observable along the geodesic 7 
and thus Hp is as well, and third, Eq. (3.28) implies 

(R F ) j 0j0 = 0. (3.30) 

That is, the second spatial derivatives of the CFC metric correction [Eq. (3.15)] are trace-free 
and correspond to a pure tidal field. Eq. (3.30) of course only follows from Eq. (3.29) in the 
absence of velocity shear and vorticity. More generally, our condition Eq. (3.29) implies 

a F 2 (R F ) j 0j0 = - aij(T zj + ujiju lj . (3.31) 

So, one way of interpreting our choice of ap is that it brings the Raychaudhuri equation into 
the simple form Eq. (3.31), which algebraically relates the trace part of h F 0 to the velocity 
shear and vorticity. 

Note that Eq. (3.29) only determines of up to a multiplicative constant, which can be 
trivially absorbed in a constant rescaling of the spatial coordinates x l F . 

3.4 Residual gauge freedom at 0[(x l F ) 2 ] 

Even when combined with a prescription to fix ap(Tp), the CFC construction does not 
univocally determine the CFC metric perturbations h F v . The situation is in fact analogous 
to the case of FNC (see App. D of [11]). The metric components h F 0 and h F { are completely 
fixed at 0[(x F ) 2 ] by the requirement of conformal flatness at 0\x l F \. However, a residual 

'Note that it is not strictly a local observable, since it requires a time integration. The true local observable 
is (M/dtF ■ Nevertheless, it can be reconstructed by a local observer if we follow his geodesic through time. 
Our choice neglects any initial peculiar velocities of the geodesic congruence, which decay oc a -1 due to the 
Hubble expansion. 
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gauge freedom exists for the spatial components hf. Consider the following (time-dependent) 
reparameterization of x l F : 

x l p —> x l F + -A l jki ( rp ) xfx F x l F + 0[(x^) 4 ], (3.32) 

where the coefficient tensor A l jki ( t f) is an arbitrary function of the CFC time and is fully 
symmetric with respect to its last three indices. This change of coordinates leads to the 
following change in hf: 

hfj )> hfj + x F x l F + 0[(x l F ) 3 ]. (3.33) 

This residual gauge can be used to bring hf into a desirable shape. As an example, one can 
think of the different ways to parametrize spatial hypersurfaces of constant curvature, e.g. 
conformally flat and stereographic projections. Note that A l j kl has to drop out from the 
final result for any given physical observable. We have verified that this is the case for the 
application we will consider in Sec. 7. 

4 Perturbed FLRW spacetime 

The CFC construction described in the previous section is completely general, i.e. it can be 
applied to any spacetime. However, the most obvious application is the case of a perturbed 
FLRW spacetime, given in some global coordinates by 

ds 2 = a 2 (T)[q^ + h^dx^dx" . (4.1) 

For full generality we do not make a specific gauge choice for the metric perturbation h Fl/ . 
However, in the following we assume that is small, so that it suffices to work at linear 
order in h. The perturbative expansion in h^ u should not be confused with the power 
expansion in x l F . The former expansion is valid as long as \hfu,\ <C 1, and adopted here 
merely for calculational simplicity, whereas the latter expansion is valid if \x l F \ is smaller 
than the typical scale of variation of h^, and is the true inherent regime of validity of the 
CFC frame. We leave a generalization to nonlinear order in h„„ to future work. In realistic 
cosmological settings, metric perturbations are present on all scales, so that a coarse-graining 
needs to be performed in order for the CFC to be valid over a finite region. We return to 
this issue in Sec. 5. 

Consider a free-falling observer traveling along the time-like central geodesic. Her 4- 
velocity can be parameterised as = a (1 + hoo/2, V r f ), where the 3-velocity V* is con¬ 
sidered as first-order perturbation. The corresponding tetrad is 8 

(eo) M = a -1 ^1 + ~h 0 o,Vf\ , (4.2) 

(e,)^ = a 1 ^V 0 ,j + hoj, 6 i j — —h l j + , (4-3) 

where £ijk is the three-dimensional Levi-Civita tensor. Here and in the following, latin 
indices i,j, k,... are raised and lowered with following standard practice in cosmology. In 

8 The expressions here for the tetrads, as well as other quantities V „, ■ • ■, are not restricted to the central 

geodesic, as one can think of imaginary free-falling observers following through every time-like geodesic in the 
Universe. Therefore, partial derivatives of them are meaningful. 
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the presence of vector metric perturbations, it is necessary to introduce the rotation uj 1 . Since 
the tetrad is parallel-transported along the central geodesic, V„ and cu l obey the following 
equations, 


K + UV l Q 


LO 


kf 


-d l hoo — h l q — Hh 1 o, 

(4.4) 

£ klj (dih 0 j - djhoi) • 

(4.5) 


A prime denotes derivative with respect to the conformal time. 

We now determine the CFC scale factor ap through the definition Eq. (3.29). In the 
perturbed spacetime Eq. (4.1) d is given by 


d 


3 H + — Tihoo + —h + d ■ V 0 


(4.6) 


where h = 5 13 hij. It is understood that the right hand side is computed in global coordinates 
on the central geodesic. The local Hubble parameter is then given by 


cif{tf) 


di{r) + \u(r)h m + \ h ' + \ d ' V ° 


Moreover, the “cosmic acceleration” is 


(4.7) 


1 ddH-F^jF ) 

cl 2 f (t f ) dr F 


cl 2 f (tf) drp \a F {T F ) drp 


da F {rF)^j 


1 ' H f (tf) ^ 
3 cif(tf) 


+ ^(eo )pd^d. 


(4.8) 


Even after the condition Eq. (3.29), there is a residual freedom in of, corresponding to 
a constant rescaling of —> cap (equivalent to an integration constant when integrating the 
Hubble rate). This (trivial) residual freedom can be most conveniently removed by requiring 
ap to asymptote to a when compared at fixed CFC proper time tF , 


of(tf(If)) 


1, 


as t f —y 0 . 


a(r(t = t F )) 

For a matter-dominated universe at early times, this implies 


°f(tf) v n 1 u 

/ \ ^ t Q U00,ini; 

a(r) 3 


as tf —> 0 . 


(4.9) 


(4.10) 


where hoo.ini = ^oo(Ani) is the metric perturbation evaluated on the geodesic at early times 
on superhorizon scale. Eq. (4.7) yields 


d clf{tf ) apdi of77f drp 

dr o(r) a a dr 


ap 

a 


\h' + \d-V Q 
b 3 


and a direct integration gives 


Qf('Tf) 

a(r) 


1 g H ^00,ini “1“ 




+ 3^ 


(4.11) 


(4.12) 


This defines the physical, locally inferred CFC scale factor for a general linearly perturbed 
FLRW metric. 
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Our construction of the CFC is almost identical to that of [11], with one important 
difference in the definition of ap. Here, we have constructed ap as (essentially) a local 
observable, while in [11] was chosen to be equal to the background scale factor a, i.e. 
ap(tp) = a(t = tp). This choice however does not correspond to an observable scale factor 
and Hubble rate Tip (since differs for different gauge choices in global coordinates). In 
this case, the CFC metric corrections oc R^ kl/l in Eqs. (3.15)-(3.17) are not guaranteed to 
be observable by themselves. On the other hand, by ensuring that clf is observable, these 
corrections are observable, which is a crucial advantage. Nevertheless, the choice made in 
[11] was sufficient for their conclusions, since their goal was to show the absence of long- 
wavelength corrections to the CFC frame at lowest order in derivatives (not suppressed by 
kp). On the other hand, for long-wavelength tensor perturbations, which is the subject of 
[16] and Sec. 7, ap is equal to a at linear order in the tensor amplitude. 

5 Coarse-graining and Einstein equations in the CFC frame 

Having discussed the definition of the CFC frame and metric for a perturbed FLRW space- 
time, in this section we show how this construction removes gauge artifacts and simplifies 
the calculation of nonlinear long-short mode coupling in the regime of practical interest in 
cosmology. 

The corrections h^ v to FLRW in the CFC metric are composed of second derivatives 
of the global metric perturbations a 2 /!^ multiplied by ( x l F ) 2 . Specifically, we have terms of 
order 

tfv ~ {Kvi nh '^ d ktipvi dkdih^XpXp ■ (5-1) 

The corrections h^ v are locally observable, as the CFC construction makes explicit (note that 
a physical definition of a locally measurable scale factor is crucial for this). Thus, the equiv¬ 
alence principle guarantees that for a single Fourier mode kp of an adiabatic perturbation, 
scales as k\ in the limit kp —>• 0. Note however that individual contributions to in 
some global coordinate system are not in general observable by themselves; terms of the type 
' Hh 1 listed above are an example. Higher order corrections to the CFC metric which we 
have dropped are correspondingly suppressed by higher derivatives of h /w multiplying higher 
powers of x l F (Sec. 3.2). This means that the CFC metric is only valid over comoving spatial 
scales much smaller than Xp = kj 1 . 

In our universe, metric perturbations exist down to very small scales. For the CFC 
metric to be valid over a finite region, we need to introduce a coarse-graining of the metric 
on some comoving spatial scale L = A -1 . The CFC metric is then constructed with respect 
to the coarse-grained metric perturbation which in Fourier space has contributions only 
for wavenumbers k < A. In this section, we describe the coarse-graining and derive the 
structure of Einstein equations in the CFC frame. We work throughout to linear order in the 
coarse-grained metric perturbation, but allow the small-scale perturbations to be nonlinear. 

5.1 Coarse graining 

Consider the Einstein equations written in some global coordinate system: 

Gfiu[g\ = SttGT^ w . (5.2) 

We assume everywhere that a cosmological constant is defined into (or equivalently 
Tpu) for simplicity. Specializing to the case of a perturbed FLRW metric [Eq. (4.1)], we 
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subtract the homogeneous background solution, to obtain the Einstein equation for the metric 
perturbation h^ u : 

G a W[h}=8irGT lxu , (5.3) 

where G a G) [, h ] is the Einstein tensor for the perturbation to the FLRW metric, which depends 
on a(r). T^ v now only contains the perturbation to the stress energy tensor. Since we no 
longer need the homogeneous background stress energy tensor, we have kept the same symbol 
Tfj, v for simplicity. 

We now coarse-grain the stress-energy tensor and the metric on a spatial scale A -1 . 
This means we have to introduce a slicing of spacetime. The choice of slicing influences the 
results obtained in the coarse-grained CFC calculation, but it does not affect the structure 
of the equations, which is our main focus here. For simplicity, we will assume that the 
coarse-graining is performed on a constant-coordinate-time slice. Hence, we write 

hfiv = v + h^ u ] = T^ v + T^ u , (5-4) 

where the “s” components are defined as the difference between the full quantity and its 
coarse-grained version, i.e. through Eq. (5.4). Let us further separate the Einstein tensor 
into linear and nonlinear pieces in the metric perturbation, G a ^[h] = G( lin )[h] + G^[h], 
where is the part of the Einstein tensor that commutes with the coarse graining. 

From coarse-graining Eq. (5.3), we obtain 

G^[h A ] + G^ A [h] = 8; tGT% , (5.5) 

which can be turned into an implicit equation for /i A , 

G' 1 i" l ['! A ] + G<"„»[/> A ] = 8irGT A „ + {cgV] - Gj,»J> A [fc]} . (5.6) 

The term in curly brackets on the r.h.s. is an effective source in the Einstein equation for h A . 
It encapsulates the backreaction of small scales, which are contained in the full h = h A + h s , 
on large scales. In the following, we will neglect this term, since we only work to linear order 
in h A , and we are not interested in the backreaction effect of small scales on large scales 
(which in general is a difficult problem). 

Let us now subtract Eq. (5.6) from Eq. (5.3), dropping the additional source, to obtain 
an equation for h s \ 

G^ [h s ] + Gff [h] - G$> [h A ] = 8t rG , (5.7) 

We are interested in the mode-coupling between long-wavelength modes h A and short scales 
h s . Suppose we have a solution h to the full Einstein equations Eq. (5.3). The small-scale 
part of the solution h s can be expanded as power series in h A : 

h s = h s0 + h sA + h sA2 + • • • , (5.8) 

where h sA " is proportional to (h A ) n . Thus, h s0 is the small-scale metric perturbation that 
remains in the limit of vanishing long-wavelength metric perturbations h A —?• 0. In keeping 
with the linear treatment of h A , we will truncate the series at h sA in the following. Note that 
T s is defined through Eq. (5.4). The stress-energy is of course a nonlinear function of the 
fluid variables and the metric. For example, a component such as pg^u in T^ u is given by 

{pg^Y = pg^u - (pg^) A ■ (5.9) 
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Similar relations hold for say pu^u v . Thus, T^ u is to be expanded in a power series of the 
type Eq. (5.8) up to the same order as h s . At zeroth order, the Einstein equation becomes, 
not surprisingly, 

[h s0 ] + G$ [h s0 ] = 8t rG T$ , (5.10) 

while at linear order in b A we obtain the equation for h sA : 

G^[h sA ] + G^'[h s0 , h A + h sA ] = 8v tGT sA , (5.11) 


where the “response” 




OM-,,, M = 

1,1 dh„„ ' 


+ a(a«v) “ 


<«<"»[/.] n a I 

+ 9 (W»V) W "' 


(5.12) 


This operator is obviously linear in its second argument and corresponds to the linear term 
in the Taylor series of G^ around h (we have dropped the tensor indices for clarity). As 
desired, Eq. (5.11) is linear in b A and h sA , but in general it is nonlinear in h s0 . This is the 
equation one is solving when computing the impact of a linear large-scale perturbation h A 
on small-scale perturbations in General Relativity. 


5.2 Transforming to CFC frame 

Eq. (5.11) was derived in some global coordinate system. Our goal now is to obtain the 
analogous equation in the CFC frame. Specifically, we construct the CFC frame for the 
coarse-grained metric. Of course, the full Einstein equations Eq. (5.3) transform covariantly. 
However, we have performed a coarse-graining and “A — s” split which breaks covariance. 
Let us thus look at the transformation more carefully. The transformation to CFC 

x —> xp{x ), (5.13) 


is by assumption a “long-wavelength” one, since it is purely determined by the coarse-grained 
metric. That is, we can write 


d x n 

dxi 


= 5^ + A^ Q [h A ], 


where A is linear in /i A . Any symmetric tensor transforms as 


(5.14) 



dx a dx@ 
dx F dx v F a,i ° J 


Gfiv + 2A“ At G ai/ ) + 0(A ~), 


(5.15) 


where G F denotes the tensor in the CFC frame, and we can neglect the 0(A 2 ) part in keeping 
with our linear treatment in h A . In particular, applying this to the long-wavelength metric, 
we have by construction 

dx^ c)x^ 

~dx^ ^ 3 = 0,2 + 2A mi/ ] = a 2 F \rj^ u + h A i F ] , (5.16) 

F F 

where h A lf oc x l F x^ F in the vicinity of xp = 0. It is then easy to see that each term in 
Eq. (5.7) and Eq. (5.11) transforms in the same way, so that the equation for h sA in CFC 
reads 


G (lin)F [h^ A ] + C (nl)F, [h s0 , hf} + C (nl)F/ [h s0 , h%] = 8t tGT sA ’ F , (5.17) 

where we have suppressed the tensor indices for simplicity. Here, the tensors Cr( lin ) F and 

G (nl)F, 

are contracted for the CFC metric, in particular they involve the local scale factor 
ap in the background. Thus, the Einstein equation for hpS in the CFC frame has the same 
structure as the one written in global coordinates Eq. (5.11). 
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5.3 CFC Einstein equations 

Let us now investigate Eq. (5.17), focusing on the l.h.s. for the time being. In the following, 
we will only deal with the CFC metric. Therefore, for the sake of clarity, we will omit the “F” 
label (referring to CFC) everywhere, except in hp, which differs crucially from h A in global 
coordinates. Suppose further that we are aiming to solve the Einstein equation on the central 
geodesic xp = 0. We will return to this point below, but note that we are free to construct the 
CFC around any timelike geodesic. We will study two regimes in which the analysis simplifies 
in different ways. First, we consider the case of a second order treatment in cosmological 
perturbations, corresponding to the leading order contribution to the mode-coupling term 
h sA . Second, we will consider the general nonlinear case in small-scale perturbations, but 
restrict to the case where the small-scale fluctuations are far inside the Hubble horizon. Since 
in the standard cosmology perturbations are small outside the horizon and become nonlinear 
only well inside it, this completely encompasses the regime of structure formation. 

5.3.1 Second order in perturbations 

Restricting to second order in perturbations, the second term on the l.h.s. of Eq. (5.17) is 
higher order, as h sA is already second order. Let us then consider the third term on the l.h.s.. 
Since we need at least two spatial derivatives on hp to obtain a nonzero contribution on the 
geodesic, only a small specific subset of terms contribute to G'( nl ) , [/i s0 , hp] in Eq. (5.12). We 
can write Eq. (5.17) as 

CFC, 2nd order: Gffl [/i sA ] + O{h s 0 didjh^)^ = 8nGT/ A . (5.18) 

Note that since the Einstein tensor is linear in second derivatives of the metric, h s0 has to 
enter without derivatives here. In order to solve the Einstein equation at second order in 
CFC for Dtp = 0, we thus only have to work out a very specific type of terms, i.e. those with 
two perturbations and two spatial derivatives acting on one of them. 

One caveat to this statement is that this assumes we do not act with spatial derivative 
operators on the Einstein equations. In practice, we might want to do that to extract the 
longitudinal part of the ij component of the Einstein equations. In that case, we need 
to include the contributions of 0[x l ] and 0[(x 1 ) 2 ] to the CFC Einstein equations, while any 
higher powers of x l are consistently dropped in the CFC expansion (Sec. 3.2). This diminishes 
the calculational simplicity of the CFC Einstein equations when ks is superhorizon. 

5.3.2 Subhorizon limit 

Let us now consider the subhorizon limit for the small-scale perturbations h s , i.e. 

As always, we allow the long-wavelength perturbation h A to be arbitrarily long. In the 
subhorizon limit, we approximate 

didjh s ~ 1, dih s ~ e s , h s , h s ', ■ ■ ■ ~ e 2 , (5.19) 

and neglect terms that are of order e 2 and higher. Intuitively, this corresponds to and 
expansion in e s = 71/ks , where ks is the wavenumber of small-scale perturbations. Then, 
the leading terms in the Einstein equations are order 1 (such as, for example, the d 2 (j) term 
in the 00 component of the Einstein equation in Newtonian gauge). Following the same 
arguments as above, i.e. using that the Einstein tensor is linear in second derivatives, both 
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G( nl ) F/ [/( s0 . lt' sA ] and G^ nl ) F/ [/i s0 , hp\ appearing in Eq. (5.17) are at least order e 2 . Thus, for 
ks H and dropping terms of order s 2 , Eq. (5.17) reduces to 

CFC, subhorizon: G^ Fn,SH ^ [/i sA ] = 87rGT*^ , (5.20) 

where is the linear Einstein tensor around the local FLRW background of^f) in 

the subhorizon limit. In other words, the long-wavelength modes do not appear on the l.h.s. 
of the CFC frame Einstein equations in the subhorizon limit. Note that this is not the case 
in global coordinates, where terms of the form h A ddh s0 are in general present. 

Let us briefly consider the r.h.s. of the Einstein equations, i.e. the long-short coupling 
piece T^ A of the stress-energy tensor perturbation. The metric enters only without deriva¬ 
tives, so that /i A does not appear. Thus, only the fluid perturbations, and h s0 , fi sA can 
appear on the r.h.s. of the Einstein equations (moreover, /i sA enters only through the linear 
Taylor term (dT[h s0 ]/dh^h^). 

To summarize: in the subhorizon limit of the CFC frame, the long-wavelength modes 
do not appear in the Einstein equation, but rather enter in the small-scale dynamics only 
through the fluid equations. If one is interested in the leading mode-coupling contribution 
(second order in perturbation theory), the linear Einstein equations for h sA are thus sufficient. 

5.3.3 Examples 

Finally, let us briefly discuss two applications. First, consider the coupling of a long- 
wavelength tensor mode 7 ^ to small-scale scalar modes. We now allow for ks to be com¬ 
parable to T-L. At second order, it is easy to see that 77 - does not enter the 00, 0 i and the 
trace part of the ij component of the Einstein equations for h sA [Eq. (5.18)] at all for scalar 
perturbations in conformal-Newtonian gauge. This is because the only 3-tensor available to 
contract with 7 ^ is 5ij, which gives zero. The tensor mode does enter the trace-free part of 
the ij component of the Einstein equations, and thus contributes to an effective anisotropic 
stress at second order. In order to extract the longitudinal part of these equations, we need 
to act with the operator didj/d 2 on the ij Einstein equation. As discussed above, this means 
we need to consider 0(x l ) and 0([x 1 ] 2 ) in the ij Einstein equation. We discuss these issues 
in detail in Sec. 7 and App. D through App. F. 

As a second application, consider a long-wavelength scalar perturbation. This is the 
focus of the companion paper [29], and we only give a brief outlook here. Let us work 
in conformal-Newtonian gauge and denote the long-wavelength potential with <f>, while the 
short-scale potential is denoted as f>. As we will see in [29], all isotropic effects are ab¬ 
sorbed in the locally measurable scale factor ap and curvature Kp (correspondingly, the 
long-wavelength contribution to the stress energy is absorbed in the effective local “back- 
gound” stress energy). What remains at our disposal is exclusively the purely anisotropic 
(i.e. trace-free) tidal tensor (didj — (1/3 )dijd 2 )& from long-wavelength perturbations. Note 
that this tidal field has the exact same form as the well-known subhorizon limit, even though 
here we have not restricted the long-wavelength perturbation to be subhorizon. 

This means that for an isotropic long-wavelength perturbation, at second order the linear 
Einstein equation for the small-scale modes is exact on all scales. This can be seen as proof 
of the “separate universe” picture. For an anisotropic perturbation, the long-wavelength 
perturbation can only enter the trace-free part of the ij Einstein equation as in the tensor 
case, since as before in conformal-Newtonian gauge 5ij is the only quantity to contract with. 
However, this contribution is suppressed on subhorizon scales ks 3> H in keeping with 
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Sec. 5.3.2, so that for subhorizon small-scale modes, the Einstein equations do not contain 
the long-wavelength modes explicitly at all. 

5.4 CFC fluid equations 

In addition to the Einstein equation discussed above, we also need to consider the fluid 
equations, 

= 0 , (5.21) 

which along with the equation of state close the system. In the same manner as described 
above, we can derive the fluid equation that describes the evolution of the long-short mode 
coupling contribution T sA . Naively, hp enters into this equation with only one derivative, 
leading to a vanishing contribution on the central geodesic in CFC. However, spatial deriva¬ 
tives of the velocity diV ] ~ an appear in the fluid equations. Thus, in order to close the 
hierarchy, we need to take one further derivative. Since the velocity itself depends on <9 /i a , 
the equation for v l involves two spatial derivatives of the metric. In this way, hp enters 
in the small-scale dynamics in CFC through the fluid equations. Another way to say this is 
that we need to know the velocity at order xf away from the central geodesic in order to 
derive the evolution of the density on the central geodesic. This will become apparent in the 
explicit calculation of Sec. 7. 

6 Relating CFC to observations 

We now discuss how calculations of observables in the CFC frame can be related to mea¬ 
surements made by a distant observer. This is a crucial step in connecting calculations to 
actual measurements (for example, the statistics of galaxy clustering). Since the focus of this 
paper is the CFC frame, and the detailed treatment of this step depends somewhat on the 
observable considered, we restrict ourselves to an outline of the general procedure here, and 
leave the application to specific observables for future work. 

The mapping from CFC to observations made by a distant observer are often referred 
to as “projection effects”. They correspond to measurable quantities from the point of view 
of the distant observer (after all the events being measured are on the past lightcone of the 
observer). However, they are not locally measurable for observers comoving on the central 
geodesic of the given CFC frame. Thus, the CFC frame separates locally measurable and 
projection contributions to any given observable in an unambiguous way. We consider this 
a key advantage of the CFC formalism. The projection effects further separate into two 
contributions, which we consider in turn. 

6.1 CFC as Lagrangian coordinates and mapping to Eulerian frame 

As we have discussed, the CFC can be constructed for any given geodesic, which defines the 
spatial origin of the coordinate system at all times. Now consider an ensemble of geodesics 
of the cosmic fluid. We can label each geodesic with a continuous three-dimensional index 9 
q. This q is then simply the Lagrangian coordinate of the fluid element moving along the 
geodesic labeled by q. One possible choice for q is the spatial coordinate in some global 
coordinate system at a fixed proper time tp = const (e.g. at early times). This shows that 

9 We are considering some finite spatial region here, not necessarily the entire Universe. Also, recall that we 
are constructing the CFC for a coarse-grained metric, so that crossing of geodesics can be avoided by choosing 
an appropriately large smoothing scale. 
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the CFC frame is related to the standard Lagrangian coordinates often used in large-scale 
structure calculations [25, 26]. Since q usually refers to the position at an early time, where 
the initial conditions of cosmological perturbations are set, these coordinates are especially 
useful when connecting observables to the statistics of the initial conditions. 

However, observations of large-scale structure are made in a Eulerian frame, which is 
a frame where coordinates refer to fixed spacetime points rather than fluid elements. Thus, 
we need to transform observables given in the CFC frame to a Eulerian coordinate system. 
Suppose we are measuring some observable, e.g. the density of matter or tracers, in a finite 
region of space. For example, we could be interested in the two-point functions of galaxies 
or matter in that region. Let us assume that this region is defined on a surface of constant 
proper time tp. We can construct the CFC frame for each geodesic threading this region, and 
calculate e.g. the matter density field in that frame (Fig. 2). However, each geodesic carries a 
slightly different CFC coordinate system, while an observer charting the entire region would 
refer the density measured at each point to a common “Eulerian” coordinate system. Of 
course these coordinates are arbitrary, but in this context it is most natural to choose the 
CFC constructed for the center of the region as our common coordinate system (see the 
following subsection). Let us call this frame {x_f} and the geodesic for which this frame is 
constructed G. We will also refer to it as “central CFC frame”. There are now two options 
to deal with observables on different nearby geodesics: 

1 . Calculate the observable for each CFC frame {.Ti?/} throughout the region separately, 
and then transform the observable to the {x_f} frame. 

2 . Calculate the observable directly in the {xi?} frame, but away from the geodesic G, i.e. 
for x l F 7 ^ 0 . 

Either choice can be adopted and has to lead to the same result, since it is merely a coordinate 
choice. In the first option, one has to calculate a coordinate transformation between two 
nearby CFC frames. The second option requires keeping at least one order of x F higher in 
the Einstein and fluid equations, since we need to solve them away from the central geodesic 
of the CFC. In general, many new terms beyond those discussed in Sec. 5 will appear in the 
CFC frame away from the geodesic, i.e. at linear order in x l F . Some of these new terms 
will involve higher derivatives of the long-wavelength metric perturbation /i A , which should 
be dropped as they are higher order in the CFC expansion. Other terms will correspond 
to residual gauge freedom of the CFC construction at higher order in xp. Nevertheless, 
the equivalence of the two approaches mentioned above shows that, when performing a 
consistent expansion in derivatives of h A , the only physical effect encoded in the Einstein and 
fluid equations for x l F 0 is the effect of the relative displacement of neighboring geodesics 
(geodesic deviation) at that order in derivatives. In practice, the first option is much simpler 
to pursue (and more physically transparent), since it corresponds to merely performing a 
coordinate transformation, which turns out to be very simple. We thus follow that approach 
here. 

We want to derive the transformation from the CFC constructed for a geodesic G 
which we shall call {xj?/}, to the frame {x_f} around G. We correspondingly denote the local 
scale factors as a F > and a.p, respectively. First, let us derive the equation for the origin of the 
{xf'} frame in the { x _ f } coordinate system, which we will denote as s. Since the origins of 
both frames follow geodesics, this is at leading order given by the geodesic deviation equation. 
Note that the displacement s is not a strictly locally measurable quantity; however it is in 
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Figure 2: Illustration of the mapping between nearby CFC frames {xf}> { x F'} (tetrads 
indicated as arrows) constructed for the geodesics G, G', respectively (thick solid and dashed). 
The initial separation at tf = 0 of the two geodesics in the {xf} frame is given by q (dotted 
arrow), which at late times is modified to q + s due to the local geodesic deviation [Eq. (6.7)]. 


principle observable as long as two observers on G and G' can communicate over the distance 
s, unlike the “absolute” displacement between global Eulerian and Lagrangian coordinates 
which is not an observable. We know that in the absence of CFC metric corrections, which is 
equivalent to Rf = 0 on the central geodesic, the separation vector remains constant since it 
is a comoving separation. The effect we are interested in is thus due to a nonzero Rf, and we 
use the geodesic deviation equation for the conformal metric gj^ u . As long as we can assume 
the scale factor to be the same cif(tf ) everywhere, this gives the desired result. Below we 
will show that this is justified. One subtlety arises because the trajectories we are considering 
are geodesic with respect to the full metric, which implies they are not geodesic with respect 
to the conformal metric (although this would be true if we were considering null geodesics). 
By using the transformation of the Christoffel symbols under the conformal transformation 
[Eq. (3.5)], and rescaling the affine parameter of the geodesic (proper time) tp to A defined 
by dX = Op 2 dtp, the geodesic equation for the full metric can be rewritten as the following 
equation describing non-geodesic motion with respect to the conformal metric: 


d 2 x M ^ dx a dx^ 
~dP + ap ~dA^X 


9a/3 


dx a dx& 
~cLX~dA 


g^d-y In op, 


( 6 . 1 ) 


the deviation from geodesic motion being described by the term on the r.h.s. Along the 
central geodesic, W = 0, and this equation simplifies significantly: 

= clf'Hf ^o , ( 6 - 2 ) 

dX 


where we have defined the tangent vector to the geodesic = dx^/dX = a 2 F (e o) M = 
ap(l, 0,0,0). The geodesic deviation equation can easily be generalized to the case where 
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the curve is not a geodesic (see e.g. the derivation in [15]). Denoting the displacement as s a , 
it is given by 

+ (0.3) 

This formidable-looking equation again simplifies greatly in the CFC frame, since all the 
covariant derivatives reduce to simple derivatives with respect to r: 

OjF (aFS a, y = a 2 F (RF) a 00fl + <5“ 0 ,s 0 {opHp)' , or 

a F \a F s a ')'- 6% s° = (R F ) a 00 , s'* • (6.4) 

Before solving this equation, let us first look at the time component. This involves the 
component (Rf)°oo^ = — (Rf)oooh (recall that the indices are lowered with the conformal 
metric, g F v = rj^p). By construction of the CFC, this vanishes at all times: since the 
Christoffel symbols vanish on the central geodesic, Rp is only composed of second derivative 
terms of the metric, which are non-vanishing only for two spatial derivatives. We only have 
one spatial index at our disposal in (Rf)oooh however, so that this component has to vanish 
on G. The temporal displacement simply obeys 

a F \a F s 0, y — s° + Hpj = 0 . (6.5) 

This is a second-order homogeneous equation for s°. It encodes the fact that rp = const sur¬ 
faces are not proper time surfaces in CFC. Note however that there is no explicit dependence 
on the long wavelength mode so that this temporal displacement is not of interest for our 
purposes. Similarly, s° does not influence the spatial displacement s* since (RfY ooo vanishes 
for the same reason. We can thus consistently set s° to zero. 

Integrating the spatial component of Eq. (6.4) yields at first order in Rp 

r p dr' r' - 

s\r F )= dT"a F (T")(R F y 00jq i 

Jo a F(T ) Jo 

i f-Tp j / rr' 

= - ——r / dr"ap{r") ( didjh F 0 ) q 3 , (6.6) 

1 JO a F\F ) Jo 

where q is the initial displacement. In other words, q is the Lagrangian coordinate assigned 
to the geodesic G' in the region considered, where the central geodesic G corresponds to 
q = 0. In the second line we have expressed the result in terms of the 00 component of the 
CFC frame metric, which involves the same components of the Riemann tensor [Eq. (3.15)]. 

Note that an expansion in Rp is an expansion in the leading corrections to the CFC 
metric, in which we retain only linear order terms throughout. Hence, a term of order ( Rp ) 2 
is consistently neglected throughout this paper. Now, we note that going away from G, the 
departure of the scale factor a pi from a F is, by construction of the CFC, at least of order 
(s*) 2 . The tetrad (e F ) can be obtained from the tetrad (e F ) for the {rp} frame by parallel 
transport along the x l F coordinate lines with respect to the conformal metric g^ u = 

(this is because the CFC requirement is to preserve the form of the conformal metric). Again, 
since vanish along G, the difference between (e F ) and (e^) is at least order (s*) 2 , where 
as we have seen s l is of order R F . Thus, at the order in derivatives that we work in, it is 
entirely sufficient to only consider the effect of the centroid shift s l . 


- 23 




This leads to a trivial transformation of quantities calculated on a given geodesic within 
the region to the fixed {x F } frame (Fig. 2). For example, let pF(q) denote the density (or 
really any other local observable) calculated in the local CFC frame for the geodesic G' 
which is labeled by the Lagrangian coordinate q defined above. Then, at a given time Tp, 
this geodesic is located at a spatial coordinate Xf in the central CFC frame given by 

-I rTp J / rr' 

x l F = q l + s*(q, t f ) = q l + - / ——- / dT"a F (r") (didjh £,) q> . (6.7) 

* Jo a F\J ) Jo 

We can thus relate the density in the local CFC frame PF(q)( r F’) to the corresponding result 
in the central CFC frame pf(x.f,tf) through 


Pf(xf,Tf) = Pf(c${tf) 


= [l - S*(xf, Tp) di] PF(q) {t F ) 


L x - . . , , , (6.8) 

q=x F -s(x F ,r F ) q=x F 

where the derivative acts on the Lagrangian coordinate q. The first equality uses the fact 

that we only need to consider a spatial shift in q, while the second equality expands to 

linear order in the shift. This transformation now clearly looks like a transformation from 

Lagrangian to quasi-local Eulerian coordinates. 


6.2 Photon propagation effects 


We now have an expression for some observable, say pf(x.f,tf ) given in the CFC frame 
of the central geodesic G of the region considered. In order to complete the connection to 
observations, we need to relate this CFC frame to the observables of a distant observer, 
specifically the photon momentum (frequency or redshift, and angular direction) detected at 
the observer’s location. The mapping from a local Fermi Normal Coordinate (FNC) frame, 
constructed for a spatial region much smaller than the horizon, to photon observables is 
well understood; specifically, it was given as “standard ruler perturbations” in Ref. [13] (see 
also [27] for a review) at linear order in metric perturbations for a general perturbed FRW 
metric. Thus, as long as the small-scale perturbations within the CFC patch are far inside 
the horizon (which is the most realistic case), we can make use of these results to relate to 
observations, by first transforming CFC to FNC and then mapping FNC to observations. 

The CFC are a generalization of the FNC, and it is straightforward to derive the map¬ 
ping from a given CFC frame to FNC. Essentially, it only amounts to going from comoving 
coordinates (tf,xf) to physical coordinates ( t F ,XF ) defined through dtp = apdrp + 0(x F ), 
and x l F = cif x l F + 0(x F ). The line element in CFC is of the form 

ds 2 = a F (rp) [r]^ + h^ u ] dx^dx F , (6.9) 


where the tidal deviation //C, is 0[(x l F ) 2 ] (note that h F also includes the contribution from 
spatial curvature). 

We now first neglect the quadratic terms h^ w and construct the usual FNC (following 
Sec. 2) from the homogeneous flat FLRW metric, Eq. (2.4) with K = 0, which yields 


ds 2 = - 1 - 


1 S 2 \ ,,2 


of dtp 


f F dt F + 


r 1 %f t - -■ r -2 \ 

Oij T X 2~ \ F-iXF,j OijXpJ 

o CLp 


dx F dx 2 F , (6.10) 


where r F = 5ijX F x F . We then restore the tidal deviation h^ v , 


ds 2 = - 1 + 


+ 


1 dUp 
clf dt F 


f p ~\ —2 ~hi 


CL i 


00 


dtp T — dtpdxp 

CLp 


1 U 2 1 

dij + -— {xF,iXF,j — dijVpJj H g -hfj 

o dp i dp 


dx l pdx J , 




( 6 . 11 ) 
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Factors of 1 /a 2 F have been added, because x F in FNC refers to the physical separation while 
in CFC it refers to the “comoving” separation. A few comments on this result are in order. 

First, Eq. (6.11) is exactly the FNC metric one would construct directly from a per¬ 
turbed FLRW spacetime, after inserting the expression for the local expansion rate and 
acceleration, i.e. Hf/clf and (1 /a 2 F )dHF/<1 tf, in terms of background quantities and per¬ 
turbations, using Eq. (4.7) and Eq. (4.8). Thus, in addition to extending the regime of 
validity to beyond the Hubble horizon, the CFC provides an unambiguous separation of the 
quadratic metric corrections in the standard FNC construction into terms that describe an 
effective, locally measurable “background” expansion (which to zeroth order in h^ u is the 
actual background), and terms which cannot be absorbed into the background cosmology 
and therefore describe a departure from a local FLRW universe. In order to completely 
prove this statement, we also need to unambiguously define the effective curvature Kp in the 
CFC frame, which is not made explicit in Eq. (6.11). We will show this in the companion 
paper [29]. We just mention here that the curvature can be uniquely defined in terms of 
the spatial CFC metric perturbation hfj. This separation of terms into local Friedmann and 
non-Friedmann terms can help to clarify the physical impact of long-wavelength modes (e.g., 
[ 28 ])- 


Second, if our goal is to connect to observations, then there is no need to insert the 
expressions for clf^f and work out Eq. (6.11). As mentioned above, the mapping from 
the local FNC of the emitting source to the observer’s coordinates of observed redshift and 
arrival direction of the detected photons is already known. Furthermore, for subhorizon 
perturbations (at leading order in e s , Sec. 5), this mapping reduces to the well-known lensing 
and redshift-space distortion effects. We then only need to take into account that in CFC, 
comoving coordinates are defined with respect to clf rather than a. Thus, spatial coordinates 
need to be rescaled by 


x G = 


qf(Lf) 

a(t F ) 


x F . 


( 6 . 12 ) 


where tp is the proper time at emission (note that the mapping from observed redshift at 
observation to proper time at emission is part of what we call the mapping from FNC to 
observations). In addition, conformal time intervals need to be transformed according to 


cItg = a 1 dtp = —drp . (6.13) 

a 

No such transformation is needed for physical time intervals dt G = dtp of course. Once the 
observables are transformed through Eqs. (6.12)-(6.13), one can then apply the mapping to 
observations derived in [13]. 

For clarity, let us summarize the three components of the projection effects necessary 
to map local CFC-frame observables to measurements by a distant observer: 


1. Transform from Lagrangian coordinates to a common, “central” CFC frame {x^?} 
[Eq. (6.8)]. 

2. Rescale comoving coordinates from local scale factor clf to the observer’s background 
a [Eqs. (6.12)—(6.13)] . 

3. Perform the mapping from local FNC to observations by adding the ruler perturbations 
of [13]. 
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This three-step process may seem fairly complicated. However, it is important to em¬ 
phasize that this is a completely general scheme which works for any observable, i.e. we 
do not have to re-derive these steps for each new observable considered. Second, each of 
these three steps has a clear physical interpretation on its own and is very easy to perform 
in practice. 

7 Application: long-wavelength gravitational waves 

As specific application of the CFC frame, we consider a FLRW universe perturbed by large- 
scale gravitational waves, or tensor metric perturbations (see [14] for previous related work). 
This amounts to setting /loo = hoi = 0 in Eq. (4.1), and requiring that hij = 77 is traceless 
ij = 0 , and divergence-free d lr )ij = 0 (note that in the notation of the previous section, 
this is h A ). The CFC frame can be used to describe the local, measurable impact of the 
tensor modes over a spatial extent much smaller than the typical wavelength of the tensor 
modes; specifically, if the long-wavelength mode is a Fourier mode Ul and the small-scale 
perturbations considered have a wavenumber kg, then the terms neglected in the CFC frame 
calculation are suppressed by an additional factor of kg/ks■ In large-scale structure appli¬ 
cations, this is an excellent approximation for a primordial gravitational wave background, 
which is the case we will consider here. 

Since tensor modes are traceless, it is easy to see that the local scale factor ap [Eq. (4.12)] 
as well as the Hubble rate are the global ones; in particular, the global coordinate time t 
is equal to the proper time for comoving observers. Let us consider small-scale scalar per¬ 
turbations (j), -0, which we describe in conformal-Newtonian gauge (this is h s in the notation 
of Sec. 5). We will assume a matter-dominated universe. By construction, the CFC is not 
restricted to subhorizon scales; we thus do not require ks 'H. For simplicity we will re¬ 
strict to the leading order effect on small scales, and thus insert the linear solution for the 
small-scale fluctuations when calculating the coupling between long wavelength tensor and 
short scalar modes. 

In the following, we shall drop the subscript p from all quantities for simplicity, since 
we will deal exclusively with CFC frame quantities. In that frame, the metric becomes 

ds 2 = a 2 (r) [(—1 — 2(f) + Loo) dr 2 + 2ho; drdx 1 + ((1 — 2ip) dij + h ij) dx l dx^] , (7.1) 

where the local influence of the long-wavelength 77 is encoded in the CFC tidal terms h 

hoo = Ij (Wki + Yki) x k x\ (7.2) 

hoi = ^ {da'ki ~ dklii) x k x\ (7.3) 

hij = (-7 ij tiki ~ Yki tiij + Yki tiij + ikj tin) X k x l 

+ ^ {dkdi jij + didj 7 ki ~ didj 7 k i ~ d t di 7 k j) x k x l . (7.4) 

Here, the coefficients multiplying x k x l are to be computed along the central geodesic. Since 
we assume negligible first-order anisotropic stress, 7 ij satisfies the equation of motion 7 " + 
21/77 ~ < 9 2 7 ^ = 0. 

In the remainder of the section we proceed as follows. We begin with the conceptually 
simpler subhorizon case H in Sec. 7.1 (recall that throughout kg can be arbitrarily 
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small). We then describe how the results of the CFC calculation can be related to observations 
in Sec. 7.2. Finally, we drop the approximation kg Ti and solve the full second-order 
Einstein-fluid system in Sec. 7.3. 

7.1 Subhorizon short modes 

We first restrict the short modes to be in the subhorizon regime A ;g H, but otherwise 
leave the tensor wavelength fully general. As discussed in Sec. 5, it is sufficient to retain 
the Einstein equations at linear order in the small-scale modes. These assume their familiar 
form in conformal-Newtonian gauge, where we can moreover identify the two potentials 
ip = (j) on subhorizon scales. The Einstein-fluid system for the short-scale fluctuations is 
most conveniently described by the subhorizon Poisson equation, continuity equation and 
Euler equation, 


d 2 <p = -U 2 5, 

(7.5) 

5' + di [(1 + <5y] = 0 , 

(7.6) 

v[ + Uvi + (v J dj) Vi = -dip + ^ ( 7 '(j + U'y'ij) x 3 , 

(7.7) 


Note that while gravity part is linear, we retain the full nonlinearity in fluid motion. The 
CFC tidal terms do not enter the first two as we restrict to x l = 0. However, the Euler 
equation must be kept at 0\x l \ to account for a velocity shear djV 1 at x l = 0, and the long- 
wavelength tensor enters at that order by supplying a tidal force (l/2)<9jhoo- We also define 
the velocity divergence 9 = diV 1 . 

In the following, we will solve the system Eqs. (7.5)-(7.7) to second order. Following 
our discussion in Sec. 5, this is not a necessary restriction: one could solve the Einstein-fluid 
system to fully nonlinear order in small-scale modes in the subhorizon regime, e.g. using 
an N-body simulation, without having to worry about nonlinear terms involving the tensor 
mode in the Einstein equations. 

Following our perturbative approach, the above equations are solved as follows. We start 
with solutions for <5, v and (p linear in scalar and tensor perturbations, at which order scalar 
and tensor perturbations are decoupled. We then compute corrections to short-wavelength 
scalar perturbations from their couplings with long-wavelength tensors. With this procedure, 
the solution, taking 5 for example, is 

<5 = h (1) + 5 (2) , (7.8) 

where <5^ is the linear solution, 6^ is the second-order correction due to mode coupling 
(denoted by a superscript sA in Sec. 5), and correspondingly for the other variables (p, v l . 

For adiabatic initial conditions, and neglecting the decaying mode for the scalar contri¬ 
bution , 10 the linear solution is given by 

<5 (1) = ^-2 + ^ 2 d ^j ‘/’ini, E (1) = - + ^7 ij x 3 , # (1) = ~^d 2 (p ini . (7.9) 

Throughout matter domination the potential remains constant ^^(x, r) = 0i n j(x). The 
linear density and velocity divergence only need to be known right on the central geodesic, 

10 This is justified since the scalar fluctuations are of much smaller scale than the tensor modes, and thus 
enter the horizon long before the tensor mode does. 
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while the velocity itself needs to be accurate to O [x l ]. It is worth highlighting that in CFC a 
contribution to v l linear in the long-wavelength tensor is induced at 0\x l \. At second order, 
Eqs. (7.5)-(7.7) read 


5 2 0 ( 2 ) _ 3 H 2 5 (2) 

5®' + e® 
0 (2)/ + H 0 (2) + S 2 c/> (2) 


= 0, 

= -di 




- (< d iV (1)i ) - (v^di dj u (1)i 


(7.10) 

(7.11) 

(7.12) 


When linear solutions are inserted on the right hand sides, we only have to keep mixed 
quadratic terms proportional to linear short-scale scalar perturbation multiplying linear long- 
wavelength tensor, as we focus on tensor-scalar coupling. An equation for 5^ can be then 
derived by eliminating 0 ^ and <9 2 0^. When combining equations, we take time derivatives, 
but additional spatial derivatives are avoided in order not to mix up different orders in x z in 
CFC. We find 

s (2) " + ns®' - ^n 2 sw = - . ( 7 . 13 ) 

Before integrating this equation, we need to specify initial conditons. This corresponds to the 
physical coupling of long-wavelength tensor modes with short-scale scalar modes when the 
latter cross the horizon during inflation. Since the tensor modes have long exited the horizon 
at that time, we do not expect the scalar modes to show any physical coupling with these 
long-wavelength tensor modes. Indeed, [11] showed that the squeezed-limit tensor-scalar- 
scalar bispectrum in single-field inflation vanishes when evaluated in the CFC frame. Thus, 
the initial condition is (5( 2 ^(r; n i) = 0 = ^^^(rini). The solution reads 


SW( T ) = --^-S N (k L ,T) 7ini didj<j>i n i = -S N (k L ,r) 7 ^ \t) 


d 2 


rkr.T 


S N (k L ,T ) = - / dx 


dTy(x) 

dx 


1 - 


k L r 


(7.14) 


Here is the initial, superhorizon amplitude of the tensor mode and 7iy (kr) = 3ji (kr) / (kr) 
is the tensor linear transfer function. In the second equality, we have used that the linear 
density constrast in the synchronous-comoving gauge 5$ is related to the Newtonian gauge 
potential in Einstein-de Sitter universe through <5^ = (2/3'H 2 )d 2 (j)i n i. 


7.2 Relating to a distant observer 

Eq. (7.14) describes modified clustering due to long-wavelength gravitational waves as a local 
observer would see at her location. A distant observer, however, would measure how density 
varies over a finite region of space. If the distant observer wants to compute the power 
spectrum of the apparent clustering, she would need the leading order correction to account 
for deviation from the CFC central geodesic (x l 7 ^ 0) in which a result like Eq. (7.14) is 
derived. Fortunately, as already discussed in Sec. 6.1, this does not require an independent 
calculation that goes to the next order in x l , but merely a spatial shift as in Eq. ( 6 . 8 ). 

Applying Eq. (7.2) to Eq. ( 6 . 6 ), we find the shift vector to be 


s* = 2 Q 3 


dr’ 

a(r') 


J dT " a ( T ") [7 ij + (! - T y ) , (7.15) 
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To leading order, we replace the Lagrangian coordinate q l with x l and apply Eq. (6.8) to find 
the quasi-local Eulerian density up to second order 


5 e = 5 {1) + ^ (1 - %) 7-ni Xi djd {1) - S N (k L ) 


a(i) 

Q2 sc 


(7.16) 


On subhorizon scales, which we have assumed throughout this section, <5^ = and we 
have complete agreement with the expression obtained in [16] [Eq. (55) there]. 

Turning to the comparison with [20], the second term in Eq. (7.16) manifests itself as 
the first term in Eq. (34) there; note that in that calculation the bispectrum in comoving 
gauge [1] was used as initial condition in the second order calculation, leading to the factor 
of (1 — Ty) of this term. The third term is equal to the second term in Eq. (34) of [20] 11 . 
The third term in the latter is due to the mapping from local FNC to observations, i.e. the 
third contribution in the list of projection effects at the end of Sec. 6 (the second in that 
list being absent in this case, since clf = a). If we include this projection effect as well, in 
the notation of [20], we obtain for the anisotropic distortion of the density power spectrum 
[defined in Eqs. (31)—(32) there] 

^ +B « Ax » - . 

(7.17) 

where P$ is the linear matter power spectrum. Note that when performing a Fourier transform 
of <5 e, the explicit dependence on x l present in CFC turns into a derivative with respect to 
ks (see e.g. App. F of Ref. [20]). 

The CFC thus automatically provides a clear physical interpretation of each of these 
contributions, which makes it considerably easier to check the calculation for correctness. In¬ 
deed, even in the very restricted case of long-wavelength tensors coupling to short-wavelength 
scalars, gauge artifacts can be very difficult to control when calculating in global coordinates 
at second order. 


7.3 Short modes on all scales 

Since CFC is not restricted to subhorizon distances, we now generalize beyond the limit 
ks T~L assumed in Sec. 7.1. We derive the fully general relativistic result for the second- 
order density measured on the central geodesic. Compared to the subhorizon case, this 
does involve some additional subtleties in the CFC resulting from the longitudinal part of 
the trace-free ij component of the Einstein equations. First, while the 00, Oi, and trace 
component of the Einstein equation can be taken to be linear as in the subhorizon case, the 
trace-free component contains terms of the type d 2 ^ and ’Hyh multiplied by (f >, ?/>. Second, 
unlike the other Einstein equations, we need to act with a spatial derivative operator in order 
to extract the longitudinal part of the trace-free component. 

It is convenient to use the 00-Einstein equation, the divergence of the Ck-Einstein equa¬ 
tion (with linearised Einstein tensor), supplemented with the (nonlinear) equation of energy 


11 The sign of this term in Eq. (34) of [20] is a typo. 
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conservation, as well as the divergence of the momentum conservation, 


d 2 *i> - - 3h 2 </> - \n 2 b = o, 

(7.18) 

d 2 (/ + n<j>) - ^n 2 di [(1 + s) v*] = 0, 

(7.19) 

S'+ di [(1 + <5K] -3(1 + <S)^ = 0, 

(7.20) 

O' + 'H9 + d l [{v-’dj) + d 2 (j) = 0. 

(7.21) 


According to Sec. 5, these equations are valid at x z = 0. At first sight, these are four equations 
for four unknowns <j), i\>, 5, 0 , and should thus suffice to solve the system. However, this is not 
correct in CFG due to the breaking of translation invariance: in order to solve the system, we 
would have to take inverse spatial derivatives, which cannot be consistently done at a single 
point but requires higher orders in x l to be included as well (see App. D). One can however 
close the system without taking additional (inverse) spatial derivatives of any of the above 
equations by including the longitudinal part of the trace-free Einstein equations, as we will 
see below. 

Note that we now have to distinguish between the two scalar potentials <and i)j, since 
at second order they are no longer equal (although their difference is negligible on subhorizon 
scales ks S> 77). We therefore introduce Physically, this stems from an 

effective anisotropic stress induced at second order by linear perturbations. Following the 
perturbative procedure adopted in Sec. 7.1 [where the linear solutions still given by Eq. (7.9)], 
we obtain equations for second-order perturbations, 


S 2 V> (2) - 377^ (2)/ - 377 2 </> (2) - -77 2 <5 (2) 

-d 2 (V 2 )'+?^( 2 )) _^ 2 0 (2) 
6®' + e® - 3V> (2)/ 
0 {2)/ + ne {2) + <9 2 </> (2) 


o, 

3 

2 


V 2 d i 



-& 




+ 3 <S( 1 >^ (1)/ , 


(7.22) 

(7.23) 

(7.24) 

(7.25) 


Again, on the right, only quadratic source terms proportional to the product of linear scalar 
modes and linear tensor modes need to be kept. Without taking additional spatial derivatives, 
we can combine Eq. (7.23) and Eq. (7.25) to find an equation for d 2 ^M 2 ' along G, 

(d 2 ^ 2 )]' + m ( d 2 ^ (2) y = -n-fk +nd 2 T^'. (7.26) 

Alternatively, we can combine all equations Eqs. (7.22)-(7.25) to obtain an equation for ^ 2 ^ 
itself along G, again without taking spatial derivatives, 

V> (2) " + 3 77 ^ (2)/ = 1 < 9 2 r( 2 ) + 77 r^'. (7.27) 

3 

Comparing Eq. (7.26) and Eq. (7.27) illustrates the point made above: we cannot simply 
take the Laplacian of the latter to obtain the former, since we have neglected terms of 0[x 1 } 
and 0[(x 1 ) 2 ] in both equations. 

Let us now turn to the equation for (9 2 r( 2 ) at x l = 0, obtained from the longitudinal 
part of the trace-free ij component of the Einstein equations, that is, 

/ Q^Qi 1 ■ A 

[-ffr--/ 1 )^Gi3-^GT l3 \=W (7.28) 
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In order to consistently evaluate this equation including the derivative operator acting on G ^, 
we need to explicitly evaluate Gij to second-order (retaining terms of scalar perturbations 
multiplying tensor perturbations) and up to 0[(a7) 2 ] (since any power of x l from h^ u can be 
converted into an inverse derivative acting on the short-scale perturbation, see App. D). We 
then set x l = 0 at the very end. This is technically somewhat involved, but can be done as 
detailed in App. E. We finally obtain on G 


0 2 T (2) 


x ( 7 ij + + <9 2 7 ij) 


d i d j 

d 2 




(7.29) 


Apart from 9 2 T^ 2 ), we also see that T^ 2 )' enters Eq. (7.27). However, unlike d 2 T^ 2 \ terms 
of r( 2 ) and r( 2 )' are negligible at the order in CFC expansion that we work in. To see 
this, note that the solution for T^ 2 ^ will involve more inverse derivatives acting on (j)^\ More 
specifically, a net negative power of derivatives act on . These terms cannot be determined 
unambiguously from the CFC metric up to 0[(a7) 2 ], as they are sensitive to 0[(x 1 ) 3 ] terms in 
the CFC metric, which we have neglected throughout. They are also affected by the residual 
gauge freedom Eq. (3.33) up to 0[(x 1 ) 2 ]. Physically, this means that these terms are beyond 
the leading local effect of long modes in terms of an expansion in ki,/ks- Thus, we neglect the 
term in Eq. (7.27). We then integrate Eqs. (7.26)-(7.27) to find 

<9V 2) = -SAT^Kni ini) + 3 Sat(/cl) ’ ( 7 -30) 

V )(2) = -Sfebini ) , (7.31) 

where 


<5 (/cl) = 5at(/cl) + 1 - %f{ki), (7.32) 

as introduced in [20]. From Eq. (7.22), we then find the density on G at second order 
5f = <5^ + 5^ to be 


5f 


«<■> + ~ 2 (^ <2) ' + ^ 1> ) + 2r<2) 


2 -\ , 2 
— 2 + , n 8 I (f )ini — 


+2 7i 


377 2 

13 I 

ini 1 Q2 ' 


ini 


3 ^ 2 57v(A:L)7ini {9idj<t> ini) + ^S N {k L ) [8 7j 

1 dS(k L ) 


77 dr 


+ 5(/cl) 



(7.33) 


We have neglected the term involving T* 2 ) in the first line because the leading result for it 
involves a net negative powers of derivatives acting on q i^ 1 ). Eq. (7.33) shows that the long- 
wavelength tensor modes couple to small-scale scalar modes only through didjcj) in conformal- 
Newtonian gauge, which is directly related to the density perturbation in synchronous gauge. 

As shown in App. F, the full general relativistic result for the density Eq. (7.33) agrees 
with [20], where the computation is carrried out in global coordinates, once a coordinate shift 
relating the two gauges is taken into account. 12 Note that when ks is of order 77, the second 

12 Note also that [20] neglected the term corresponding to T-Ld 2 Y since it is numerically suppressed by 
(fcn/fcs) 2 after the tensor equation of motion is applied. However, this term is an unambiguous prediction 
from the CFC metric up to 0[(x 1 ) 2 ] and should be kept. The results agree when this term is included. 
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order density given in Eq. (7.33) which is calculated in conformal-Newtonian gauge is not an 
observable, hence results in different gauges do not need to agree. However, as we show in 
App. F, the conformal-Newtonian gauge in CFC, and the global conformal-Newtonian gauge 
are at lowest order related by a simple coordinate shift related to the initial tensor mode 
amplitude. Thus, the agreement in the second order density is a nontrivial check of the two 
different calculations. Finally, it is worth pointing out that the case ks ~ T~L is not really 
of practical interest, since it implies that the tensor mode with k l <C ks is far outside the 
horizon, suppressing its observable effects. 

8 Summary and conclusions 

In this paper we have presented a systematic way to construct Conformal Fermi Coordinates 
(CFC). These are the coordinates of a free-falling observer that describes her neighborhood 
as an FLRW space up to corrections that grow with the square of the distance from her 
worldline, as summarized in Eq. (3.15), (3.16) and (3.17). The coordinates are defined for 
any, arbitrarily inhomogeneous and anisotropic spacetime, as long as the geodesic maintains a 
finite distance to any singularities. In cosmological applications, one typically wants to define 
CFC with respect to some long-wavelength perturbation. The coordinates are then valid 
only in a region that is much smaller than the wavelength of that perturbation; specifically, 
we have shown that higher order corrections are suppressed by powers of the conformal 
Riemann tensor and its spatial derivatives, (d m RF) n Xp rn+ ^ with m + n > 1. The key 
advantage that makes CFC very useful for cosmological computation is that they are not 
restricted to subhorizon scales, contrary to what is the case for standard Fermi Normal 
Coordinates (FNC). Indeed, the standard coordinates chosen by cosmologists to describe 
their observations are the CFC constructed for our entire observable universe. There are 
three major perks of using CFC. First, CFC precisely describe what a local observer would 
measure, eliminating the issue of gauge artifacts at nonlinear order from the start. Second, the 
computation of non-linear mode coupling is simplified considerably with respect to standard 
gauge choices employed in cosmological perturbation theory (e.g., Newtonian or comoving 
gauges) when ks 3> Ti (the case of greatest practical relevance). We have shown (Sec. 5) that 
in that case it is sufficient to retain the linear Einstein equations for the small-scale modes. 
Finally, the CFC allow for an unambiguous connection to observations via a well-defined 
coordinate transformation (from local to global scale factor and from Lagrangian to Eulerian 
coordinates), and gauge-invariant projection effects which have been calculated previously. 
Moreover, they cleanly separate different physical contributions of the long-wavelength mode 
due to genuine local gravitational effects, local geodesic deviation, and photon propagation 
(projection) effects. 

For example, Ref. [11], where CFC were first introduced, showed in a very simple cal¬ 
culation that there is no modulation of small-scale perturbations by large-scale potential 
fluctuations in single-field inflation in the attractor regime. That is, local primordial non- 
Gaussianity is completely absent. Similarly, Ref. [11] showed that there is no physical cou¬ 
pling of large-scale tensor modes to small-scale perturbations in single-field inflation. This 
corresponds to the absence of 0(1 ) and O(ki) effects of long-wavelength metric perturbations 
of wavenumber in the regime <C 1. 

The advantages of the CFC just described now also allow us to calculate the leading 
observable effect of long-wavelength modes in single-field inflation, which starts at ( k^/ks ) 2 . 
As an explicit example of this, we have computed the effect of a long-wavelength tensor mode 
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on a short-wavelength scalar perturbation (the result is given in Eq. (7.33)). For a subhorizon 
scalar perturbation, this had already been computed in [16] using CFC, and we recover the 
same results in that limit here. A calculation in global coordinates was presented in [20], 
which in the calculation of the second order density did not assume subhorizon small-scale 
modes (the connection to observations was however still done assuming ks S> 'H). Although 
the calculation for horizon-scale small-scale modes in CFC involves some subtleties, it can be 
done and we show that we find perfect agreement with the global calculation once the proper 
gauge transformation (a simple coordinate shift) is taken into account. The conclusion is that 
the quantity capturing the mode coupling effect on all scales is the second spatial derivative 
of the potential 0 in Newtonian gauge, which is directly related to the density perturbation 
in synchronous-comoving gauge. 

Getting rid of all this apparent nonlinearity of gravity might seem like a bit of a miracle. 
However there are several reasons why this works. First, by isolating local observables, we 
remove a large number of terms that are merely gauge modes. Second, we are relying on a 
separation in scale between long- and short-wavelength modes, so that the long-wavelength 
modes are almost constant over the scales relevant for the short-wavelength modes (specifi¬ 
cally, the leading tidal contribution of the long-wavelength mode is approximated as spatially 
constant). Finally, we are making use of the fact that while the Einstein equations are lin¬ 
earized in perturbations , the background is solved to fully nonlinear order. This ensures that 
all effects of a long-wavelength mode that can be captured by an effective background of(tf) 
are solved for exactly, i.e. at fully nonlinear order. This is related to what is commonly called 
the “separate universe” conjecture, which will be studied in detail in the companion paper 

[29]- 

Finally, it is important to remember the limitations of our approach. By construction, 
the CFC only capture gravitational effects of the long-wavelength mode. Thus, they are valid 
as long as non-gravitational interactions are irrelevant on the scale of the long-wavelength 
mode, i.e. as long as the long mode is outside of the sound horizon of any non-gravitational in¬ 
teractions. Hence certain computations, specifically those where the non-gravitational physics 
is nonlocal on the scale of the horizon 'Hr 1 ^ must be done in global coordinates. A specific 
example would be the effect of long-wavelength tensor modes on scalar perturbations before 
recombination. 
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A Evaluating ap and its derivatives 

In this section we give coordinate-free expressions for the first and second derivatives of ap, 
which are necessary in order to derive the explicit transformation from some coordinate frame 
to CFC. These derivatives are already fixed by the fact that of = of(tf) in CFC and by the 
form of the metric Eq. (3.1). 
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We start with the local comoving expansion rate along the central geodesic, 

U F (P) = <lhl ^ P) = ^ (V M ln a F ) p = a F (P)(e 0 f p (V^Inap) p . (A.l) 

Since the gradient of a F along the central geodesic points in the tangent direction of the 
central geodesic, we can solve this for the gradient, 

= (A.2) 

a F {r) 

This determines the tensor of shift Eq. (3.6). In order to determine the coordinate transfor¬ 
mation up to third order in x, we also need derivatives of the conformal Christoffel symbols, 

again evaluated on the central geodesic. This involves the derivative 

C^a/ 9,7 = s a V-yV^ In a F + 5% V 7 V Q In a F - g a p g^ x V 7 V A In a F 
- (g a p„ 9 M + 9ap 9^\ 7 ) V A In a F + (s& T^ + ^T^ - g^f^x) V, hi . (A.3) 

Note that the metric Eq. (3.1) implies that the second derivatives Vo-V^ln a F appearing in 
the first line are given by the corresponding result for a homogeneous FLRW metric (spatial 
curvature does not enter at this order). This is because the conformal Christoffel symbols 
quantifying the departure of the metric from FLRW vanish along the central geodesic. Writ¬ 
ing this in a coordinate-free form finally yields 


(V q V /3 In a F )p 


1 dU F (P) ( 'H f (P) 
a 2 p (P) dr F " V a F (P) 
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(eo)a,p(eo)/3,p — 
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B Coordinate transformation 


In this section we derive the transformation from some global coordinate system to CFC 
explicitly up to third order in x l F , continuing from Eq. (3.9). Higher-order coefficients aJJ are 
recursively computed using Eq. (3.4). At second and third order, we obtain respectively 
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(B .2) 


The transformation from the global coordinates into the CFC, expanded to third order in 
Xp, is then given by 


Sx“ ee x“(Q) - x“(P) = a F (P)( e ,r p x‘ F - A I ry {a F (P)]‘ (<*)% x‘ F x’ F 


-- - 2Cf5, [« r (P)]’( El )J («,)“ (e t )’ p 444 + O (4) .(B.3) 
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Transforming the metric 


We now show how the metric transformation can be calculated explicitly using the coordinate 
shift derived above. Coordinate shifts up to 0[(x^) 3 ] as in Eq. (B.3) fix the Jacobian ma¬ 
trix dx^/dx u F to 0[(x l F ) 2 }, corresponding to the leading corrections in CFC frame from the 
background FLRW metric. Here and throughout, P denotes a point on the central geodesic. 
Evaluating Ox^/Btf requires the derivatives along the central geodesic, 
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(B.4) 

(B.5) 

(B.6) 

(B.7) 


The first equation holds since (eo )p is the tangent unit vector of the central geodesic, while 
the second defines the conformal Hubble parameter of the CFC scale factor. The third 
relation is a consequence of the stipulation that (ei) p are parallel-transported along the 
central geodesic, and the fourth is obtained from the third directly through the chain rule. 
We eventually obtain 


1 dx» 

cif{P ) drp 


= (e 0 )p + a F 


— ter - r^(eo)“(ei)* 

a,p ^ 


x l F 


J p 


-a 2 F 


+o[& F n 


pfj, I _ PM per 

aF a/3 F- I 2 «/3,7 L aa L P'y 
i \3 


(eo F (e i )%(e j f P x' F a? F 

. p 
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The computation of dx^/dx^, by comparison, is more straightforward, since all quantities 
defined on the central geodesic do not depend on x l F . Hence we have 
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following directly from Eq. (B.3). Eq. (B.8) and Eq. (B.9) can be compared to Eq. (A9) 
and Eq. (A10) of Ref. [19], respectively, which provide the corresponding relations for the 
standard FNC. Instead of the Christoffel symbols the conformal version appears r^, 
in addition to of and Tip/ap • Note that ap and x l F always appear with the same power, so 
that only the combination clfx 1 f matters. 

When Eq. (3.10) is computed order-by-order in x F , we have to perform an expansion 
at several places. First, we have derived order-by-order expressions for the Jacobian matrix. 
Second, the conformal metric g^ u = a~ 2 g iiu has to be expanded about the CFC origin, 

9^(x x ) = (g FU )p + (d a guv) P dx a + ^ (d a dggp U )p 5x a 6x 13 + 0[(5x^)\ (B.10) 
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where Sx p is given by Eq. (B.3). Last but not least, we need to expand the factor (a(Q)/a(P)) 2 
as well. Since the global scale factor is only a function of the global time r, we first find the 
global time lag of a point at ( t f ,x f ) with respect to the CFC origin at (t f . 0 ), 


St( X ‘ f ) ee T (Q) - t(P) = a F (P)(e,)° P x' F - \ (f£„) p [a F (P)] 2 (e,)J( ej )J 44 + 0[(4) 3 ]. 

(B.ll) 


This induces a change in the global scale factor 

= 1 + H(P)6t(x),) + \ (W(P) + H\P)) MPp)] 2 + 0[6t 3 (x)?)), (B.12) 

where 'H(P), 'H'(P) are to be evaluated on the central geodesic. 

An important simplification arises in Eq. (B.10) when considering a linearly perturbed 
FRW universe. The conformal metric g FV (x) = (a/a F ) 2 (g Fl , + h /lL ,(x)) deviates from the flat 
metric only at first order in h^, and the zeroth-order piece is a constant and does not 
contribute to derivatives. It is therefore sufficient to include the zeroth-order shift Sx^ = 
5 1 x l F , which is purely spatial [Eq. (B.3)]. 

Conformal flatness of CFC at 0\x l F ] 

We now explicitly verify that the CFC metric is conformally flat, up to tidal corrections of 
0[(x 1 f) 2 ]. This is an independent but equivalent calculation to the one in the CFC frame de¬ 
scribed in App. C. The strategy is to check each of the three types of components of Eq. (3.10) 
- 00-conrponent, Oi-component, and ij -component. Algebra (e.g. partial derivatives) are 
carried out in global coordinates. 

For the 00-component, we expand to sufficient orders in x l F . Useful simplifications 
include a 2 (P)(g FV ) P = (g FU )p and a 2 (P)(g FU ) P (e a ) P (epYp = g a p. In particular, we need 
expressions up to linear order in x l F , 

5x p = ap(P)(e;)^Xp, St = a F (P)(e i ) 0 x l F . (B.13) 
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The last two terms cancel because d\ In a = TiSox. The first two terms also cancel because 
2 9ap^ = dpgpv T dvgpfi — dp9 P v- 
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We then check the Oi-component. The conformal connection will be involved. We 
need the following expression at leading order in x l F , 
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The first term will cancel with the fourth term by Eq. (B.15). The second, the third and the 
last term add up to zero after expressing the Christoffel symbols in terms of d\g a p. Therefore, 
the whole expression vanishes at 0[xp]. 

Finally, we check the zj-component. We find 
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The third term cancels with the last term by 5 a hr a = P5q\. The remaining three terms sum 
up to zero. In particular, since (7“^ vanishes when contracting with three spatial tetrad 
vectors, the extra pieces in the conformal connections do not contribute. 

We therefore have proved Eq. (3.1). The proof does not assume that the global metric 
g^v is parameterised in small perturbations, and therefore holds at all orders in perturbation 
fields (not to be confused with expansion in x l F ). 


C Structure of CFC metric corrections 

This appendix supplements algebraic details underlying the discussion of Sec. 3.2. That is, 
we show that the metric in CFC is of the form Eq. (3.1), and derive the explicit expression 
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for the 0{x 2 F ) corrections. We also discuss the higher order terms that appear at x F and 
higher. Throughout we will only deal with quantities in CFC, nevertheless we will keep the 
superscript F here. 

The CFC are defined by, first, that Eq. (3.11) holds for all times along the central 
geodesic. Second, we require that = 0, which is equivalent to Eq. (3.14). To show 

this, we consider again the curve Eq. (3.12). It describes a geodesic with respect to the 
conformal metric g, and it is linear in the affine parameter A. In CFC, the tangent vector is 
(0,/?*), and the conformal geodesic equation Eq. (3.4) reduces to 


(f F )l 


ftp = o, 


(C.l) 


which implies (T^)C|p = 0 because /3 l is arbitrary. To examine the remaining symbols 
(f F )(| a , we use that (e a )^ = a^ 1 5a are parallel-transported along the central geodesic. In 
CFC, ( eo) l 'V l/ (e a ) ,J = 0 reduces to 

a~ F l d o (a^O + a F 2 (r%« = 0, (C.2) 

which becomes 
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where throughout this appendix a prime denotes a derivative with respect to t f . We then 
merely have to relate the ordinary symbol to the conformal one I r (components evaluated 
in CFC) using Eq. (3.5) and Eq. (3.6). The gradient of a F is required to be tangent to the 
central geodesic, so we have In ap = (In ap)' and 


(° F )oa = (w)'c 


Combining this with the result for (F f )q 0 leads to 


(f 


F\P 
0 a 


= o. 


(C.4) 


(C.5) 


This guarantees that there are no terms linear in x l F in the CFC metric. 

To derive the 0[(x l F ) 2 ] terms, we directly apply the strategy of Ref. [12], For FNC, 
the authors deal with Christoffel symbols and Riemann tensor with respect to the physical 
metric. For CFC, all results carry through, as long as we instead deal with quantities defined 
with respect to the conformally-related metric. To be specific, we want to relate d a dp g pi> to 
Rp vpa in CFC at some point along the central geodesic. Again, since (W )“ v always vanishes 
there, we have 


Besides, we have in CFC 


<9o(f F )% 


= 0 . 
p 


(C.6) 


^(f F )“ 0 =(R F )%» o 


(C.7) 


The remaining derivatives ^(r^)^- involve three spatial indices downstairs. The derivation 
parallels the one detailed in Ref. [12], which we do not repeat here but only outline the 
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idea. We can consider the geodesic deviaton equation (with respect to g) applied to the 
tangent vector m ~ d/dx l F of the spatial geodesic (see also Sec. 6.1; note that unlike the 
central geodesic itself, the spatial CFC coordinate lines are geodesic with respect to g). The 
geodesic deviation equation has the schematic form VVn + R F n = 0, where the second 
covariant derivative reduces to the spatial derivative of m and therefore leads to a formula 
df f = R f . The final result reads 


F'\ct 


4(r% 


(R F Ti jk + (R'Yuk 


>F\a 


J P 


(C.8) 


Combining Eq. (3.11) with Eqs. (C.6)-(C.8), we obtain double spatial derivatives of g , 


didj g F o\p 
didj <?ofc|p 
d i d j9ki\ P 


~2 R F ioj 


-v (R 


H) ikj 


+ Ri 


0 j ki 


-- R 


Hkjl 


+ r f 


jlik 


(C.9) 

(C.10) 

(C.ll) 


These results imply Eqs. (3.15)-(3.17). 

We now discuss what the next, higher order contribution to the metric at 0[(xp) 3 ] looks 
like. On the central geodesic G , we examine 


didjd k g F v \ p , (C.12) 

where P denotes an arbitrary point on G which we keep fixed. Note that we do not need 
to consider time derivatives in CFC, since these can be consistently taken within the leading 
CFC frame metric. (We hereafter drop the super/subscript F as we will deal exclusively with 
quantities in CFC). We have 

Kb = (2 d {a g vp) - d v g a e) . (C. 13) 

Thus, taking two spatial derivatives and evaluating on G yields 


d i d iKp p = if*didj (2 d {a g u p) - d„gap) 


+ 0 [(dr v )(dd~g jS ), (ddrn(d~g yS )] . (C.14) 


Only the first term, involving three derivatives of the metric, is nonzero on G. In order to 
extract the triple spatial derivative terms of the metric (as explained above, we already know 
all the other ones), we need to consider didjT As shown in [12], since I vanishes 

on G we have 

U P Ok 


— — T-t ,J 
P ~ n fco j 


(C. 15) 


We then obtain 


f) r) T^ 1 
Ok 


= -diR 1 


kOj 


(C. 16) 


since the correction from the (T) 2 term in R vanishes even after taking a derivative. The 
spatial derivatives of r(\ can be obtained by taking another derivative with respect to the 
tangent T along the spatial geodesic of the geodesic deviation equation. Again neglecting 
terms of order T dT which vanish on G, this yields simply the ordinary CFC spatial derivative 
of the geodesic deviation equation written in CFC, Eq. (64c) of [12]. Schematically, 




= d, 


K Hj+ p erm - 


(C.17) 
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To conclude, the 0[(x l F ) 3 } corrections to the CFC frame conformal metric scale as 

(d[R^ juk J p x' F a? F x k F , (C.18) 

i.e. they are suppressed by the spatial derivative in CFC of R F multiplied by xp. At order 
(x F ) 4 , we similarly have terms of order ddR F , but we also in general obtain terms of order 
(. R f ) 2 . Thus, we are performing an expansion in powers of R F and its spatial derivatives, 
with appropriate powers of xp. 


D Inverse spatial derivative in CFC 


In this Appendix we discuss how inverse spatial derivatives are to be performed in CFC, 
which break translation invariance. We show that it amounts to the following. One formally 
writes d = ds + dg and then Taylor expands in powers of dg/ds, subject to the rule that 
<9g’s only act on short-wavelength perturbations and dg's only on x l F . The algebra is then 
straightforward because no negative power of dg exists. Since we consistently drop any 
0([x *] 3 ) and higher corrections to the CFC metric, this expansion is truncated at d\. 

In CFC, at 0[(x l ) n ], a generic term on which one desires to act with (inverse-)spatial 
derivative looks like 


f(x l , t)F(t) x n x 12 ■ ■ ■ x ln , 


(D.l) 


where f(x l ,r) only depends on the short modes, and F(r) only depends on the long modes. 
Note that F(r) only depends on time by definition, since it is evaluated on the central 
geodesic. We have suppressed the spatial indices of / and F for clarity. To simplify our 
notation further, we suppress the time dependence hereafter. 

Eq. (D.l) does not have a well-defined Fourier transform, since it is not bounded at 
infinity. In order to treat this, we let F(r) have a smooth, very large-scale shape in real 
space, corresponding to very narrow support around 0 in k space. Thus, in fc-space we 
multiply F(r) with 

■ (°- 2 ) 


where C is a normalization constant and A is the cutoff in Fourier space. Apart from factors 
of order one, A is nothing else than the coarse-graining scale introduced in Sec. 5. If we take 
A —y 0, this factor becomes (2ir) 3 5 F \k), i.e. the Fourier transform of unity as desired. If 
we use ~ to represent the Fourier transform of a field in real space, the Fourier transform of 
Eq. (D.l) is a convolution 


f(x l ,T)F(r) x ll x 12 



d 3 kg C_ ~ - -id -id 
(2ir) 3 A 3 ^ dk?ldk?l 


—id 


exp 



,(D.3) 


where kg = k — kg. The number of —id/dkg operators acting on the exponential counts the 
order in x l . 

We now use d~ 2 as an example to demonstrate how inverse, and therefore non-local, 
spatial derivatives can be computed for a power expansion in x l . It amounts to power 
expanding the operator in Fourier space 


-1 
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1 - 2 gr k i + 
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K s 




+ 4 


ks,iks,j 

ikiy 


kik[ + 0[(k L y 
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so that ki, only appears in the numerator. Powers of ks then trivially translate into (inverse- 
spatial derivatives of the short mode /. 

For each power of ks, we can commute it with —id/dk^s until it appears nearest to the 
exponential; but then this term vanishes in the limit A —> 0. To convince ourselves, we have 


k l L exp 



iA 2 —id f k 2 T \ 
~Ydk^ exp {-A 2 )- 


(D.5) 


In real space, this seemingly raises the power of x % by one; however, because of the two extra 
powers of A, it approaches zero once the limit A —> 0 is taken in the end. Therefore, the 


surviving terms are those generated from the commutation relation 


ki, d/dk J L 


5 l j, 


which in real space translates into [x l , dj ] = —S l j, identical to the action of <9j on powers of 
x l . The net result of our procedure boils down to writing d l = d l s + d l L , Taylor-expanding 
in dL/ds, and following the rule that dg ’s only act on /, and d^’s act on powers of x\ For 
instance, we have for the non-local operator <9*<9 J /d 2 (which is used in our calculation for this 
work) 


d l di 
d 2 


*s4 

d 2 s 
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dr 


’S,l 


0 d 2 s ) 2 


+ 4- 


■ 9sd§d s ,kd s ,i 


C d 2 s ) 3 


d k L d l L 


+ o 



(D.6) 


The power expansion in d^/ds truncates for local operators, but is in general an infinite 
series for non-local operators. Therefore, non-local operators in principle mix up all orders 
of x l in CFC. In practice, the series is truncated at finite order. The reason is that each 
additional power of <9^ comes with an inverse ds, and is therefore suppressed by one more 
power of ki/ks- In our computation, we drop terms that are proportional to a net negative 
power of ds acting on the short-scale potential. These necessarily come with three powers of 
ds and thus require the CFC metric beyond quadratic order. Therefore, only a finite number 
of terms in the power expansion need to be considered. 


E Second-order anisotropic stress in CFC 


We start with the second order part of the spatial-spatial Einstein equation G l j 


8 t rGTtj, 


2t/> (2) " + n (4V> (2)/ + 2^ (2)/ - d 2 (V (2) 


+ d i d j (V (2) - <^> (2) ) 


+0[(ddcj>U)hF„] + <D[(d{dhp) + 0[^ ( ddhp] = 8 t tGvMvV. 


(E.l) 


On the left hand side, we first of all collect second-order short-scale perturbations. In addi¬ 
tion, from an expansion of G l j at second order, we have a bunch of quadratic terms involving 
linear short modes (or its derivatives) multiplied by the CFC tidal metric h^ u (given in 
Eqs. (7.2)-(7.4)) or their derivatives, which we do not explicitly spell out here. Since the 
CFC metric correction h^ v is truncated at 0[(x 1 ) 2 }, these terms are at most 0[(x 1 ) 2 ]. On 
the right hand side, we have a term quadratic in v \*from T l j. We only have to insert the 
solution for up to 0[x l ] (Eq. (7.9)), because the higher orders are generated by CFC 
metric correction h^ v beyond the quadratic order, which is subleading in ks/ks- 
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The longitudinal scalar part of this tensorial equation is obtained by acting upon it with 
the non-local operator 


d 2 



(E.2) 


as defined in CFC by the Taylor expansion Eq. (D. 6 ), and then evaluating at x l = 0. Accord¬ 
ing to the general discussion in App. D, it is sufficient to expand to 0[(<9 l) 2 ] in Eq. (D. 6 ). 

With the prescription of App. D, straightforward but tedious algebra can be carried 
out. However, even at x l = 0, more terms can be dropped. Indeed, we learn from App. D 
that, beyond leading order in Ol/Ds, eliminating one power of x l by one power of ds occurs 
at the expense of introducing one more power of ds in the denominator, which brings about 
an additional suppression by a factor ks/ks- Simple power counting suggests that to retain 
leading terms in ks/ks, we can drop all terms having a net negative power of ds acting on 
the short mode. 

The end result on the central geodesic reads 



(E.3) 


The left hand side comes from G l j while the right hand side comes from T l j. This leads to 
the solution for the second-order “effective” anisotropic stress Eq. (7.29). 


F Comparison with previous results 


In Ref. [20], the density perturbation to second order was computed in global coordinates, 

[T y (k L ) + S(k L )] {didj<t> in i) 


< 5 G = (-2 + ^ 2 U ini - 2 


3R 2 


+^S N (k L ) 


m 2 

didj 


d 2 


^ini ) + 27j 


didj 

~W' 


1 dS(k L ) 
% dr 


-f Slyklf) 


• (F.l) 


Only the coefficient of the second term ~ didjf i n i does not agree exactly with the CFC 
result Eq. (7.33). The discrepancy is due to the fact that although results in both frames 
depend on the same (and physical) initial potential profile 0im, d 2 (/) i n i are in general different 
scalar functions for different parameterisation of the spatial coordinates. This is because 
d 2 = 5^ didj is not coordinate invariant. 

If the spatial origin of both frames is chosen to coincide, then to leading order, on the 
initial spatial slice, the global parameterisation can be converted into the CFC parameteri¬ 
sation (proper distance) via a shift 


i i] 

X ~ 2^ini x 3 ■ 


(F.2) 


If we label quantities in CFC by F and those in the global frame by G , then on the initial 
slice we have (note that 7 R refers to the value at the origin and therefore should be regarded 
as a constant) 
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(F.3) 
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Taking the spatial derivative twice, and then evaluating at x l = 0, we find 


(d 2 (j) ini) G = (d 2 (j)i n iY - 7 y i (didj<f>i D i) F ■ 


(F.4) 


At the order of perturbation we pursue, this only applies to the first term (linear solution) 
of Eq. (F.l) and yields an additional term at second order, so that 


Sg 


2 ,,\ 2 
~2 + — , 0 8 ) </>i n i — 


3H 2 


+^S N (k L ) (d 2 ^ 


3U 2 

didj 


[Ty(k L ) - 1 + S(k L )] 7 A (didj(j) i n i) 


Q2 <t> ini J + 27 ini ^ g 2 y ml 

Given Eq. (7.32), this is in complete agreement with Eq. (7.33). 


didj 


1 dS(k L ) 
H dr 


+ S(k,L ) 


.(F.5) 


References 

[1] J. M. Maldacena, JHEP 0305, 013 (2003) [astro-ph/0210603]. 

[2] P. Creminelli and M. Zaldarriaga, JCAP 0410, 006 (2004) [astro-ph/0407059]. 

[3] A. Kehagias, H. Perrier and A. Riotto, Mod. Phys. Lett. A 29, no. 29, 1450152 (2014) 
[arXiv:1311.5524 [astro-ph.CO]]. A. Kehagias, J. Norea, H. Perrier and A. Riotto, Nucl. Phys. 
B 883, 83 (2014) [arXiv:1311.0786 [astro-ph.CO]]. P. Creminelli, J. Norea, M. Simonovic and 
F. Vernizzi, JCAP 1312, 025 (2013) [arXiv: 1309.3557 [astro-ph.CO]]. M. Peloso and 

M. Pietroni, JCAP 1305, 031 (2013) [arXiv:1302.0223 [astro-ph.CO]]. A. Kehagias and 
A. Riotto, Nucl. Phys. B 873, 514 (2013) [arXiv:1302.0130 [astro-ph.CO]]. P. Creminelli, 

J. Gleyzes, M. Simonovic and F. Vernizzi, JCAP 1402, 051 (2014) [arXiv: 1311.0290 
[astro-ph.CO]]. 

[4] C. Cheung, A. L. Fitzpatrick, J. Kaplan and L. Senatore, JCAP 0802, 021 (2008) 
[arXiv:0709.0295 [hep-th]]. 

[5] P. Creminelli, J. Norena and M. Simonovic, JCAP 1207, 052 (2012) [arXiv: 1203.4595 [hep-th]]. 

[6] K. Hinterbichler, L. Hui and J. Khoury, JCAP 1208, 017 (2012) [arXiv:1203.6351 [hep-th]]. 

[7] L. Senatore and M. Zaldarriaga, JCAP 1208, 001 (2012) [arXiv:1203.6884 [astro-ph.CO]]. 

[8] K. Hinterbichler, L. Hui and J. Khoury, JCAP 1401, 039 (2014) [arXiv: 1304.5527 [hep-th]]. 

[9] W. D. Goldberger, L. Hui and A. Nicolis, Phys. Rev. D 87, no. 10, 103520 (2013) 
[arXiv:1303.1193 [hep-th]]. 

[10] X. Chen, H. Firouzjahi, M. H. Namjoo and M. Sasaki, Europhys. Lett. 102, 59001 (2013) 
[arXiv:1301.5699 [hep-th]]. 

[11] E. Pajer, F. Schmidt and M. Zaldarriaga, Phys. Rev. D 88, no. 8, 083502 (2013) 
[arXiv:1305.0824 [astro-ph.CO]]. 

[12] F. K. Manasse and C. W. Misner, J. Math. Phys. 4, 735 (1963). 

[13] F. Schmidt and D. Jeong, Phys. Rev. D 86, 083527 (2012) [arXiv: 1204.3625 [astro-ph.CO]]. 

[14] K. W. Masui and U. L. Pen, Phys. Rev. Lett. 105, 161302 (2010) [arXiv:1006.4181 
[astro-ph.CO]]. S. Dodelson, E. Rozo and A. Stebbins, Phys. Rev. Lett. 91. 021301 (2003) 
[astro-ph/0301177]. 

[15] S. Carroll, Spacetime and Geometry: An Introduction to General Relativity , Addison-Wesley, 
2003. 

[16] F. Schmidt, E. Pajer and M. Zaldarriaga, Phys. Rev. D 89, 083507 (2014) [arXiv:1312.5616 
[astro-ph.CO]]. 


-43 









[17] T. Baldauf, U. Seljak, L. Senatore and M. Zaldarriaga, f JCAP 1110, 031 (2011) 
[arXiv:1106.5507 [astro-ph.CO]]. 

[18] C. Wagner, F. Schmidt, C. T. Chiang and E. Komatsu, arXiv:1409.6294 [astro-ph.CO]. 

[19] F. Schmidt and D. Jeong, Phys. Rev. D 86, 083513 (2012) [arXiv:1205.1514 [astro-ph.CO]]. 

[20] L. Dai, D. Jeong and M. Kamionkowski, Phys. Rev. D 88, no. 4, 043507 (2013) 
[arXiv:1306.3985 [astro-ph.CO]]. 

[21] R. Bean and O. Dore, Phys. Rev. D 69, 083503 (2004) [astro-ph/0307100]. 

[22] P. Creminelli, G. D’Amico, J. Norena, L. Senatore and F. Vernizzi, JCAP 1003, 027 (2010) 
[arXiv: 0911.2701 [astro-ph .CO]]. 

[23] D. H. Lyth, K. A. Malik and M. Sasaki, JCAP 0505, 004 (2005) [astro-ph/0411220], 

[24] C. Pitrou, X. Roy and O. Umeh, Class. Quant. Grav. 30, 165002 (2013) [arXiv:1302.6174 
[astro-ph.CO]], http: //www.xact. es/xPand/index.html. 

[25] F. R. Bouchet, R. Juszkiewicz, S. Colombi and R. Pellat, Astrophys. J. 394, L5 (1992). 

[26] T. Buchert, A. L. Melott and A. G. Weiss, Astron. Astrophys. 288, 349 (1994) 

[astro-ph /9309056]. 

[27] D. Jeong and F. Schmidt, arXiv:1407.7979 [astro-ph.CO]. 

[28] B. D. Sherwin and M. Zaldarriaga, Phys. Rev. D 85, 103523 (2012) [arXiv:1202.3998 
[astro-ph.CO]]. 

[29] L. Dai, E. Pajer and F. Schmidt, work in preparation. 


- 44 - 


